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Abstract 

A fundamental property of QCD is the presence of the chiral 
anomaly, which is the primary component of the 7r° — > 77 decay ampli- 
tude. Based on this anomaly and its small (~ 4.5%) chiral correction, 
a firm prediction of the 7r° lifetime can be used as a test of QCD at 
confinement scale energies. The interesting experimental and theo- 
retical histories of the 7r° meson are reviewed, from discovery to the 
present era. Experimental results are in agreement with the theoreti- 
cal prediction, within the current (~ 3%) experimental error; however, 
they are not yet sufficiently precise to test the chiral corrected result, 
which is a firm QCD prediction and is known to ~ 1% uncertainty. 
At this level there exist experimental inconsistencies, which require 
attention. Possible future work to improve the present precision is 
suggested. 
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1 Introduction 

An important feature of Quantum Chromodynamics (QCD)[U [2], the ac- 
cepted theory of the strong interactions, is the existence of the chiral anomaly [3j 
H]. An anomaly is said to occur when a symmetry of the classical Lagrangian 
is not a symmetry of the full quantum mechanical theory. The it — > 77 decay 
is perhaps the best example of a process that proceeds primarily via the chiral 
anomaly. Consequently, an accurate theoretical prediction for the ir° — y 77 
decay rate can be made with no free parameters. Precision measurements 
of this quantity thus serve as a stringent test of the validity of QCD itself. 
This represents a particularly interesting test, since at low energies QCD is 
very hard to solve because the quarks and gluons interact so strongly and 
predictions have in general to be made by the use of numerical techniques. 
The history of the ir° lifetime then is not only of fundamental importance, 
but also represents an interesting story of the ingenuity of experimental and 
theoretical physicists. 

The existence of the tt meson was first postulated by Yukawa in 1935 in 
order to explain the short range and large magnitude of the nucleon-nucleon 
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interaction [5]. Initially, the newly discovered \x meson was thought to be 
Yukawa's particle, but the muon turned out only to have only weak and 
electromagnetic interactions. The charged ir meson was finally discovered 
in a 1947 cosmic ray experiment [6]. This was followed in 1950 by a series 
of experiments that found the n° meson[71 El [101 [11], and its primary 
decay mode into two gamma rays. This last feature is closely connected with 
the chiral symmetry of QCD[T21 EE2], which makes the ir mesons the lightest 
hadrons[14j. 

During the 1950's it was discovered that the pion family was an isotriplet 
with spin and parity = - JE The pseudoscalar nature of these lightest 
hadrons, the pions [T4"]. experimentally confirmed the conjecture by Nambu [T2"] 
that the underlying chiral symmetry of nature is spontaneously broken. In 
modern terms, the QCD Lagrangian is chiral symmetric in the limit where 
the light quark masses are zero[13]. If this symmetry were to be manifested 
in the usual Wigner-Weyl fashion, each quantum state, such as the proton, 
would have an opposite-parity partner particle. Since this is not the case, 
Nambu realized that the axial symmetry is instead realized via the appear- 
ance of massless pseudoscalar mesons (now called Nambu-Goldtone Bosons) 
so that, e.g., the opposite parity partner of the proton is a multiparticle 
state containing the proton and massless "pion". This conjecture was put 
on a stronger theoretical basis by Goldstone[16]. Of course, in the real world 
pions are the lightest hadrons but have nonvanishing mass due to the explicit 
breaking of chiral symmetry, since the masses of the up and down quarks are 
small, but non-zero[T7] [181 EH]- The modern picture of the pions then is that 
they are Nambu- Goldstone Bosons in addition to being Yukawa's mesons, 
the source of the long range part of the nucleon-nucleon interaction. They 
play this role by having relatively weak interactions with nucleons in the 
s-wave (vanishing in the chiral limit when the masses of the light quarks — > 
0) but strong interactions in the p-wave channel. 

Electromagnetic effects make the 7r ± 4.6 MeV heavier than the neutral 
pion 7T°. This means that the only possible decay of the ir° is either the two 
gamma mode or the relatively small (~ 1.2%) 7e + e~ Dalitz decay mode[T9]. 
This decay, like the two photon decay of positronium, requires that the two 
photons are El and Ml, in order to carry away the negative parity of the J w = 
0~ state[15j. This means that their electric field vectors are orthogonal and 

1 For a brief history of the experiments leading to the measurement of J 7 * = CP for the 
pion family, and the connection with the two photon decay of positronium see Perkins [15] 
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has been experimentally demonstrated in the double Dalitz ir° — > e + e e + e 
decay [20]. 

The results for the ir° lifetime (and decay width) are shown in FigJTJ 
There have been different experimental methods utilized to measure the ir° 
lifetime. The first is the direct technique, in which the mean distance that 
the 7T° meson travels before it decays is observed. Within the first year of 
its discovery it was realized that t(tt°), the ir° lifetime, was < 5 ■ 10 -14 sec 
by measuring the decay distance in low energy reactions. The early series 
of direct decay distance measurements was concluded in 1963 with the first 
definitive measurement which utilized a 16 GeV proton beam at CERN with 
the result that r(7r°) = (1.05 ± 0.18) • 10~ 16 sec[2TJ. The second experimen- 
tal procedure utilizes the Primakoff [22] effect in which an incident photon 
interacts with the Coulomb field of a nucleus to produce the tt° meson. A 
measurement of the cross section combined with detailed balance yields the 
value of t(tt°). Measurements using this technique were carried out from 1965 
through 1974 (see Sec. H]). In 1985 a greatly improved version of the CERN 
direct experiment was performed |23j. The third method, published in 1988, 
involves measurement of the cross section for the purely electromagnetic two 
photon e + e~ — > 77 — > ir° process[21]. More recently, a measurement was 
performed of radiative pion decay — n + — > e + z/7 — by the n — (3 group [25]. 
Using isospin invariance, the weak polar-vector form factor contributing to 
this decay channel is related by a simple isospin rotation to the amplitude 
for 7T° — > 77, and in this way one additional experimental number for the 
7T° lifetime has been obtained (see Sec J3.5l for further discussion). These 
measurements complete the information on which the Particle Data Book 
(PDB)average is based [T4"| [26] and the results of these experiments is shown 
in Fig{TJ along with the newly performed Primakoff measurement |27j . With 
the exception of one major outlier [28], these results are in reasonable agree- 
ment with each other. At a more precise level, looking towards a test of 
the theoretical predictions [211 [301 EI] a t the 1% level, there exist differences 
between the two most accurate measurements [23l [27] . The PDB average is 
t(7t ) = (0.84 ± 0.04) • 10~ 16 sec). As will be discussed below, however, we 
believe that this error is understated by a significant factor, perhaps as much 
as a factor of two(the averaging procedure is discussed in Sec 14. 51) . 

Before the axial anomaly was understood, a standard way to calculate 
the 7r° — > 77 amplitude was to utilize the partially conserved axial-vector 
current (PCAC) condition, which relates the pion field to the divergence of 
the axial current via d^J^ = F^m^^ where J5 is the axial- vector current, 
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Figure 1: 7r° — >• 77 decay width in eV (left scale) and t(tt°), the 
mean 7r° lifetime in units of 10~ 16 sec. (right scale). The experimen- 
tal results with errors and publication dates are: 1) particle data book 
average^; 2,3,4) Primakoff experiments( 1970- 1974) [321 ESJ |33]; 5) di- 
rect method(1985)I23]; 6) e + e~(1988)P5; 7) vr/3 experiment (2009) [25] ; 8) 
new Primakoff measurement(2011)[27]. All of these experiments with the 
exception of the last one are the basis of the particle data book aver- 
age. The lower dashed line is the predictions of the chiral anomaly [3j H] 
(r(7r° -> 77) = 7.760eV,r(vr ) = 0.838 • 10~ 16 sec). The upper solid 
line is the chiral prediction and the dotted lines show the estimated 1% 
error[29l [301 E](r(vr 77) = 8.10eV, t(tt°) = 0.80 • 10~ 16 sec). For the 
relationship between r(7r° — > 77) and r(7r°) see EqJTJ 
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(/v represents the pion field, and F w = 92.21 ± 0.03MeF[I3] is the pion 
decay constant measured in the n + — > p^v^ decay rate. However, use of 
PCAC yields r(7r°) « 10~ 13 sec, a decay rate approximately three orders of 
magnitude too slow (for details see Secj3]). Note that in this procedure the ir° 
amplitude vanishes in the chiral limit, when the masses of the up and down 
quarks and the pion are set to zero. 

It was the discovery of the chiral anomaly that solved this theoretical 
conundrum [21 H]. The existence of the anomaly requires an additional term in 
the divergence of the third component of the axial current d^J^ = F n m^4 > ^ + 
{a/ir)E ■ B where a is the fine structure constant. From this additional 
term it can be seen that the ir° — > 77 decay is via El and Ml photons as 
indicated by experiment |2Uj. Note also that this additional term survives 
in the chiral limit, and is exact therein. In fact the anomaly term is the 
dominant contribution to the ir° — > 77 decay rate, which (in the chiral limit) 
has no adjustable parameters^ H] — 

r(7r°^7 7 ) = (m„/47r) 3 (a/F n ) 2 = 7.76(W 

T = h/T tot (7r°) = 0.838 10~ 16 sec 

T tot (n°) = r(7r°^77) + r(7r ^ 7 e + e-) 

r(7r° ->■ 77)r(7r°) = ^(vr -> 77)^ (1) 

where BR(tt° — > 77) = 0.9882 [13] is the 7r° — > 77 branching ratio. However, 
Eq. [T]is exact only in the chiral limit, i.e., when the u and d quark masses 
vanish. In the real world, there exist modifications and the dominant chiral 
corrections are due to the small masses of the up and down quarks and their 
difference [T71 [T8| [T3], which mixes the isotopic spin- zero rj and 7/ mesons 
into the isotopic spin-one 7r° wave function. As discussed below, this chiral 
symmetry breaking produces a ~ 4.5% increase in r(7r° — > 77) to 8.10 eV 
(t(tt°) = 0.80 • 10~ 16 sec) with an estimated uncertainty of less than 1%[29| 
|30| [3T] and it is an important goal of modern experiments to test this firm 
QCD prediction. 

The present average experimental value for 

r(7T° -> 77) = 7.74 ± 0.37 eV, r(7r°) = (0.84 ± 0.4) • 10~ 16 sec. 

given by the Particle Data Group p3] is in reasonable agreement with this 
predicted value [29| |30| I3T] . This number primarily represents an average 
of several experiments, all of which were performed before 1985 (see SecJH 
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for a complete discussion). The quoted error of 5% is most likely too low, 
since each of the experiments appears to have understated their errors and 
also, as can be seen in Fig. (TJ from the much larger dispersion between 
the different measurements (see the discussion in Sec. 5. 4). Even at the 5% 
level, however, the precision is not sufficient for a test of such a fundamental 
quantity, and in particular for the new calculations which take the finite quark 
masses into account and are accurate at the 1% level. All of the previous 
experiments were performed with experimental equipment which by now has 
greatly improved. 

In order to begin to improve this situation a modern experiment — PrimEx — 
was performed at Jefferson Lab using the Primakoff effect technique |27[ IM| 
135] . This experiment utilized tagged photons for the first time and incorpo- 
rated many accelerator and detector improvements developed over the years. 
These included CW accelerators which provide high duty cycles and greatly 
improved angular and energy resolution for the outgoing pion and decay pho- 
tons. Such improvements enabled a significantly better measurement, with 
a 2.8% overall error as shown in FigJU and yielded a result consistent with 
the chiral prediction. Still, even with this improved Primakoff measurement 
there is still considerable room for experimental improvements. A second ex- 
periment using the Primakoff effect has also been performed by the PrimEx 
group and the data analysis is in progress [35]. 

The chiral anomaly represents quantum mechanical symmetry breaking 
by the electromagnetic field of the chiral symmetry associated with the third 
isospin component of the axial current [31 0]. The ir° decay provides the 
most sensitive test of this phenomenon of symmetry breaking due to the 
quantum fluctuations of the quark fields in the presence of a gauge field. 
Considering the fundamental nature of the subject, and the 1% accuracy 
which has been reached in the theoretical lifetime prediction, it is important 
for future experiments to aim for the same level of precision. 

An interesting aspect of the history of the n° lifetime is the degree of 
independence of experiment and theory. This is even more remarkable since 
one of the early pioneers in the discovery, properties, and early theory of the 
7T° was Jack Steinberger. In most of the experimental papers on which the 
Particle Data Book average is based there is no comparison of the experimen- 
tal results with theory. This is true even for papers which appeared after the 
seminal work on the chiral anomaly [31 0]. It is only in the past decade that 
the PrimEx experiment [2 7J was designed to test QCD via the predictions of 
the anomaly plus the chiral corrections [221 EDI EI]- The chiral predictions in 
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turn were stimulated by the prospect of the PrimEx experiment. 

With this interplay of theory and experiment in mind we below review 
both the theoretical and experimental approaches to the ir° — > 77 decay. We 
begin in section 2 by examining the 1950 discovery of the neutral pion and 
its decay into two photons together with the early lifetime measurements 
which gradually converged towards 10 -16 sec. by 1963. In section 3 we 
review the theoretical evolution which led to our current understanding of 
this process. In section 4 we examine the experiments that are used by 
the PDG in computing their average, and in section 5, we look at the new 
PrimEx experiment performed during the last few years at JLab. In section 
6 we briefly examine some related issues. Finally, in section 7 we summarize 
our findings and speculate on future improvements. 

2 Early Experimental History 

In 1947 the charged pion was discovered with photographic emulsions ex- 
posed to cosmic rays at mountain altitudes[6], and its dominant, weak, muon 
neutrino decay mode tt + — > fi + + was observed. In 1950 the neutral pion 
was observed at the 184 inch Berkeley synchrocyclotron via proton bombard- 
ment of nuclei [7J, as well as in the n^p — > ir°n and ir~p — > 771. reactions with 
stopped pions, and the dominant electromagnetic tt° — > 77 decay mode was 
detected[8j |9]. The neutral pion was also detected in cosmic rays at 70,000 
feetpX]. In the same year the tt was also photoproduced at Berkeley and 
the coincidences between the two decay photons were observed for the first 
time[TU]. By the end of 1950 then the following facts were established about 
the 7T° meson: 

• The value m(7r+) - m(7r°) = 5.42 ± 1.02MeV^\ ^, consistent with the 
presently accepted number 4.59 MeV|14j. 

• The cross sections for the 7p — > ir°p (TUIES] and 7p — > ir + n^7\ reactions 
are roughly equal, indicating that the 7r° and ir + mesons are "of the 
same type" , indicating that the 7r° meson is a pseudoscalar. 

• The soft component of cosmic rays is due to the production and decay 
of 7T° mesons. 

• A lower limit for the lifetime t(7t°) < 5 • 10~ 14 sec was established by a 
measurement of the geometric size of the decay region [TT]. 
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It is impressive that, within a year of its discovery, so much was under- 
stood about the ir°, including a stringent upper limit of < 10~ 14 sec. for the 
lifetime. This value is far shorter than electronic detection resolution time 
and was obtained by setting an upper limit on the distance between the n° 
production and decay. This upper limit [TTj utilized the best experimental 
method that was available for such short lifetimes — studying 7r° production 
and decay in emulsions — since the resolutions are of the order of the grain 
size ~ 0.5/i. As the mean decay distance is d(7r°) = 7/3cr, we find, using the 
the predicted lifetime r = 0.81 ■ 10~ 16 sec, that cr = 0.025// rs 5% of a grain 
size! With the benefit of hindsight then, it is not surprising that actual mea- 
surements (as opposed to upper limits) of the n° lifetime were much slower 
in coming. 

In 1951 Dalitz proposed the existence of the 7r° — > / ~fe + e~ decay mode 
and calculated a branching ratio of ~1.2%[l9]. in excellent agreement with 
the current experimental value of 1. 174(0. 035)%[H]. Dalitz's primary point 
was that the observation of this decay mode would possibly enable a mea- 
surement of t(tt°) since the opening angles of the electron-positron pair are 
on average larger then those due to pair production of one of the decay pho- 
tons. In a cosmic ray interaction where a high energy particle creates many 
particles including a 7r° in an interaction ("star"), this leads to a radial distri- 
bution N(r) of e + e~ pairs, N(r) ~ const + (S/d)e~ r ^ d where S is the relative 
probability to produce a e + e~ pair and d(7r°) = 7/3cr is the mean 7r° decay 
distance. 

In 1953 a measurement of N(r) in cosmic rays, as suggested by Dalitz, 
was carried out [38]. The conclusion of this work was that "the most probable 
value of t(7t°) is ~ 5- 10 -15 sec and it is very improbable that the true value is 
less than 3 ■ 10~ 15 sec or larger than 10 -14 sec." [38] • Perkins has pointed out 
that this value should be corrected downward due to the reduction in ionizing 
power when the e + e~ pair are very close together [39] and this is the direction 
needed to bring this result into better agreement with the upper limits of 
2 • 10~ 15 sec and 1 • 10~ 15 sec found in cosmic ray emulsion experiments [40] 
as well as in low energy experiments on pion charge exchange at the Chicago 
cyclotron [41]. It should also be noted that the long lifetime claimed by 
Anand[38] depends strongly on the 7r° momentum distribution in the cosmic 
rays, a quantity that is not well determined in the emulsion experiments. In 
view of these issues the determination of Anand[38] must be considered as 
only a first tentative step in the road that lay ahead. 



9 



One of the main advances in this regard came through the development 
of higher energy particle accelerators, so that the intensity and control of the 
primary beam was greatly improved over the use of cosmic rays. In 1957 an 
ingenious method was proposed to measure the ir° lifetime from stopped K + 
mesons using the two-body decay mode K + — > 7r°7r + [l2]. The kaons were 
produced in the Berkeley Bevatron and were stopped in an emulsion. The 
decay location was determined from the appearance of the ir + . However, the 
emulsion is insensitive to the gamma rays from the dominant tt° — > 77 decay. 
Therefore the pair (Dalitz) ir° — > ^e + e~ decay mode, which occurs with a 
1.2% probability [H], was utilized. For a stopping kaon the pion momentum 
is 205 MeV/c and for an assumed n° lifetime of 10~ 15 sec the mean decay 
distance is 0.46 microns. Emulsions were utilized since they have a spatial 
resolution of one grain or « 0.5 microns. The experiment indicated that the 
7T° meson decayed in a significantly shorter distance so that t(tt°) < 1 ■ 10~ 15 
sec. g2]. 

Several years later, during the period from 1960 to 1963, the first definitive 
measurements of the ir° lifetime were reported. These experiments used the 
Berkeley Bevatron and the CERN SPS cyclotron, along with emulsions with 
better spatial resolution by a factor of ~ 2, as well as having better statistics. 
The results of these early measurements are summarized in Tableland FigfJJ 

The three earliest experiments that obtained results [431 EU EE) used 
the technique previously suggested [4"2"] of using stopped kaons and observ- 
ing the K 27r decay mode. To illustrate what a tour de force such experi- 
ments were, they observed a mean decay distance of 0.088 ± 0.24yU (com- 
pared to a developed grain size of ~ 0.35//) which leads to a mean life 
r = (1.9±0.5)-10 -16 secjH]- This number is reasonably close to the presently 
accepted 7r° lifetime — 0.84 ■ 10 -16 secfH] for which the mean decay distance 
is 0.038/i. The fourth experiment at Berkeley utilized a 3.5 GeV ir~ beam 
to produce neutral pions in emulsion nuclei and then observed their Dalitz 
decay [46]. This method depended on understanding the ir° momentum spec- 
trum, and the assumption was made that this was identical to the measured 
7r + spectrum. All of these tour de force emulsion measurements obtained a 
lifetime ~ 2- 10~ 16 sec. On the one hand this is a remarkable result, coming in 
the early 1960's. Nevertheless, as can be seen from Figj2j the values are higher 
than the presently accepted number and probably result from a systematic 
bias in the technique. One possibility (mentioned above) was that pointed 
out by Perkins in connection with the first reported measurement [38] — that 
the lifetime value should be corrected downward due to the reduction in 
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ionizing power when the e + e~ pair are very close together |39j. 

The next step in performing a measurement of the decay distance of the 
7T° meson was to utilize higher energies so that there exists a large Lorentz 
boost. In 1963, using an 18 GeV internal beam at the CERN proton syn- 
chrotron, the yield of 5 GeV positrons was measured from platinum targets 
of various thicknesses. In this case the ir° mesons, produced in the nuclear 
interactions of the protons, decay into two photons, some of which are con- 
verted in the target to e + e~ pairs. The 7r° decay distance, inferred from the 
relative positron yield as a function of target thickness, was determined to 
be (1.5 ± 0.25)/i [21]. To obtain t(tt°) the 7r° momentum spectrum must be 
known and, as in ref. jlB] mentioned above, this was taken to be the same as 
the measured ir + spectrum. With this assumption < p w o >, the average 7r° 
momentum that produced 5 GeV positrons in Pt, was 7.1 GeV and a lifetime 
r(7r°) = (1.05 ± 0.18) • 10~ 16 sec was obtained[2Tj. which is much closer to 
the the present values, both experiment al[H] and the theoretical, as can be 
seen from Figj2j 

The first Primakoff measurement came was performed in Frascati in 
1965[47j. This is the measurement of the cross section for the 7 + 7* — > 7r° 
reaction where one photon is incident on the virtual photon (7* from the 
Coulomb field of a nucleus (the Primakoff effect is discussed in detail in 
Sec. 5.1). In this case the incident photons were produced in an electron 
synchrotron with an endpoint energy of 1.0 GeV incident on a Pb target. 
The results of this experiment are in very good agreement with the present 
accepted and theoretical values. 

With these last two measurements we have arrived at the beginning of 
the era on which the particle data book is based. Before examining them, 
however, it is useful to review the corresponding theoretical studies of ir° — > 
77 and their connection with QCD. 

3 7T° — >■ 77 Decay: Theory 
3.1 Early History 

Just as the experimental study of 7r° — > 77 took place over many years, the 
corresponding theoretical understanding of neutral pion decay evolved over 
several decades. The theoretical examination of the decay amplitude for the 
mode 7T° — y 77 began contemporaneous with the work on renormalization of 
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Table 1: First Measurements of the it Lifetime. 



Reference 


Reaction 


t(tt )/10~ 16 sec 


number of events 


03] 


a 


3.2±1.2 


26 


m 


a 


1.9±0.5 


76 


m 


b 


1 q+1-3 


44 


m 


a 


z -°-1.0 


61 


m 


c 


1.05 ± 0.18 




m 


d 


0.73± 0.11 





Reaction footnotes: a) K + — > tt + ti observing the Dalitz decay mode 
7T° — > 7e + e~ with stopped Kaons; b) 7r~ — > tt° charge exchange reactions 
in emulsion nuclei with a 3.5 GeV pion beam; c) direct measurement of the 
7T° decay length induced by 5 GeV/c protons incident on a Pt foil (see text for 
discussion) at CERN.; d) The earliest Primakoff measurement at Frascati. 
The first four experiments were performed at the Berkeley Bevatron, using 
emulsions as detectors. 
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quantum electrodynamics (QED). In 1948 Sin-Itiro Tomonaga sent a letter 
to J. Robert Oppenheimer, who was then director of the Institute for Ad- 
vanced Study in Princeton, in which he described some of the work that he 
and his group had been doing in the area of QED. This letter was subse- 
quently published in as a letter to the editor in the Physical Review, and in 
it Tomonaga described his successful work in dealing with divergences involv- 
ing the electron mass and charge [4"8"] . However, he was still having difficulty 
dealing with the renormalization of the photon propagator in that the pho- 
ton self-energy diagram shown in Figure [3^ was not only divergent but also 
violated gauge invariance, leading to a nonzero value for the photon mass. 
In order to study such issues further, two of his associates — -H. Fukuda and 
Y. Miyamoto — undertook a calculation of the process 7r° — > 77, which in- 
volves the pseudoscalar meson-vector current-vector current (PVV) triangle 
diagram shown in Figure |3b[l9]. They also examined the axial current- 
vector current- vector current (AVV) triangle diagram in Figure |3b, which is 
a three-point function connecting the axial vector current to two photons. 
Such triangle diagrams are linearly divergent and this problem was dealt with 
by the use of a Pauli-Villars regulator [50J, [5 1J. The calculation raised inter- 
esting problems in that the PVV amplitude was found to be gauge invariant, 
while the AVV amplitude was not. 

A parallel calculation was undertaken by Jack Steinberger at Princeton 
(then a theorist), who was aware of the Fukuda-Miyamoto work and, also 
using the Pauli-Villars method, obtained similar results[S2]. Defining 

•£^77 A^yjTT F^iv (2) 

where 

Fplv = d^A v — d u A^ 
is the electromagnetic field tensor and 

is the dual tensor, Steinberger determined, using a proton loop, that 

A _ e 2 g wNN 
lbir^m N 

where g n n n is the strong pseudoscalar nNN coupling constant. Using the 
Goldberger-Treiman relation [53] 

m N g A 

9-kNN = ~B — ( 4 ) 
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(a) (b) (c) 



Figure 3: Diagrams considered by early workers. 

where g& ~ 1.27 is the neutron axial decay amplitude and F n ~ 92.2 MeV is 
the pion decay constant, Eq. [3] can be rewritten as 

_ e 2 g A 

*~ n lGvr 2 ^ 1 ' 

which is remarkably similar to the value 

e 2 

A anom _ /c\ 

-77 - i 6lT 2 Fn w 

predicted by the chiral anomaly, as shown below. The corresponding decay 
rate predicted by Steinberger — 

T. 77 = lAn £ m " = 9 2 A r~ = A x 7.76eV (7) 

is about 60% larger than the prediction of the chiral anomaly (and the PDB 
experimental value). 

Julian Schwinger also visited these problems in 1951 [33]. He showed how 
to handle the issues with the photon self-energy, and confirmed that there 
were difficulties with the triangle diagrams, but the source of the problem 
remained mysterious. 

3.2 The Adler-Bell-Jackiw Anomaly 

The resolution of this problem and its connection with symmetry was not 
understood until the late 1960's, when at CERN John Bell and Roman Jackiw 
examined the problem of 7r° decay within the a model, which is known to 
obey what we now call chiral symmetry[3]. At the time, the chiral symmetry 
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was manifested in the validity of PCAC, which asserts that the divergence 
of the axial current 

J^{x)=${x) l -r allJ rt^{x) (8) 

can be used as an interpolating field for the pion — 

d^J^(x) = F w m 2 Mx) (9) 

By this logic the divergence of the axial triangle diagram should be related 
to the pion decay triangle diagram, but this was clearly not the case since 
the PVV amplitude is gauge invariant while its AVV analog is not. By a 
careful analysis within the a model, Bell and Jackiw were able to show that 
the origin of the problem is the breaking of chiral symmetry in quantizing 
the theory That is, the chiral symmetry is valid classically but is destroyed 
via quantization — a situation which is called anomalous symmetry breaking 
or simply the anomaly. (At the same time as Bell and Jackiw were resolving 
the problem, Steve Adler at the Institute for Advanced Study came to the 
same conclusion in his study of spinor field theory jl], and for this reason the 
phenomenon is often called the ABJ anomaly.) 

The basic reason underlying this behavior is that, because quantum field 
theory involves an infinite number of degrees of freedom, the short distance 
properties of the theory do not coincide with what is suggested by naive 
manipulations. That one must be very careful in this region is suggested by 
the feature that in a spinor field theory obeys the anticommutation relation 

{ip a (t,x),4(t,y}} = 5 3 {x-y)5 ab (10) 

so that one must deal very carefully with a current such as JgJx), which 
involves both the field and its canonical conjugate defined at the same point 
and can be handled in a variety of ways, but all such methods lead to identical 
results. 

Since the result of all techniques is the same, we shall detail here only 
the most intuitive of these procedures — perturbation theory — which was the 
method employed by the early investigators. Before examining this calcula- 
tion, however, we review the simple symmetry aspects that one might expect. 
We consider a simple massless spinor field carrying charge e coupled to the 
electromagnetic field, for which the Lagrangian density is 

£ = # fi-e^-^F^ (11) 
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We note that this Lagrangian is invariant under a global phase transformation 
of the spinor field 

t/j -> exp(z/3)^ = $ (12) 
which, by Noether's theorem, leads to a conserved polar-vector current— 

J 11 = $1ij& with 8^.1^ = 0. 

Alternatively, the Lagrangian of Eq. [TT] is also invariant under a global axial 
phase transformation 

ij> -> exp(rn 5 )^ = i/j' a (13) 
which, by Noether's theorem, leads to a conserved axial-vector current 

hii = 4>1ilJ5^ with d^J^ = 
Consider now the three point AVV amplitude, which we designate by 

W?i>fc) = j dVye- 1 ^- 1 ^ < 0\T(j; m (x)Jr(y)J 5j (0))\0 > 

(14) 

where it is understood that q\ — q\ = 0. Current conservation for the vector 
and axial-vector currents then yields the conditions 

QiT^j(qu <?2) = q2 T ^j(qu 92) = (<?i + ?2) 7 ^ 7 (gi, q 2 ) = (15) 

The requirement that all three conditions in Eq. [T5] be satisfied then leads 
to the vanishing of the 7r° — > 77 decay amplitude, a result which is called 
the Sutherland- Veltman theorem [55], [56]. This result can be seen by writing 
the most general form for T^(qi,q 2 ) which satisfies the strictures of Bose 
symmetry, parity conservation, and gauge invariance — 



T ilvl (q 1 ,q 2 ) = exaa/3 Pj9 ^ A\<hG\{$> ) + [f ^qiv - 9° ^lJ q?q 2 9 ~,G 2 (p 



( (g\qiu - g\qiM<& - \p 2 9\g\{qi ~ q2f)g\G 3 (p 2 



(16) 



where p = qi + q 2 is the momentum carried by the axial current. Imposing 
the condition for axial current conservation yields the constraint 

= P 7 Vrfa, q 2 ) = e^q^p'iG^p 2 ) + G 3 (p 2 )) (17) 



17 



Defining the off-shell 7r° — > 77 amplitude as 

<77k°>=^V(?i»?2) ( 18 ) 

where 

A flv (q ll q 2 ) = A(p 2 )e liVal3 q<?ql (19) 
we have, using the LSZ reduction and Eq. [9] 

A(P 2 ) = ^"fW ifr 2 ) + G 3 (P 2 )) (20) 

whereby A(0) =0, which is the content of the Sutherland-Veltman theorem, 
and asserts that in the chiral symmetric limit, where m 2 = 0, the 7r° — > 77 
decay amplitude vanishes. Of course, in the real world m 2 7^ and we must 
extrapolate from the chiral limit. However, this scenario would suggest a 
decay amplitude of the size 

A ^ ~T&F„ x f <21) 

.2 / A 2 



where A x ~ 47ri^ ~ 1 GeV is the chiral scale|57].|58j. Here the factor m^/A^. 
represents the feature that this amplitude is two chiral orders higher than the 
vanishing lowest order term, the factor e 2 /47r is needed because we have a 
two photon amplitude with a loop diagram, and the "extra" 4ttF 7T is required 
for dimensional purposes. In any case, Eq. [21] would lead to a 7r° lifetime 

r„-o^ 77 = l/r^o^^ ~ 10~ 13 sec, (22) 

three orders of magnitude slower than that found experimentally. 

In order to study this phenomenon further we shall examine the 7r° decay 
process in perturbation theory, wherein the three point function is described 
by the Feynman diagram in Figure [4b 

T^(qi, q 2 ) = £V 7 (<?1, qi) + U v ^(q 2 , qi) (23) 

where 

rr r \ -e 2 K F f d A s m / 1 1 1 \ 

U ^ ^ = I (2^ Tr ( TVi^l^i^ ) (24) 
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Figure 4: Perturbation theory diagrams for PVV and AVV processes. Here 
the indices /i, v are Lorentz indices of vector currents while 7 is the Lorentz 
index of the axial- vector current. 
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Note that U fll>1 (qi, q-z) arises from colored quark loops and includes a factor 



K F = N c Y,Q 2 q r 3q = 3 



u.d 



(25) 



This the proportionality of the decay amplitude to N c has led many authors 
to assert that the agreement between experimental and theoretical values of 
the 7T° — > 77 decay rates offers "proof that N c = 3. However, Baer and 
Wiese have noted that anomaly cancelation requires, in a two-flavor picture 
with N c colors, that the u, d quark charges must have the values 



Q u = -[— 

^ u 2 [ N r 



1] and Q 
so that the factor Kp has the value 

H) 2 



4 [K N r 



1 



1. 1 



if 



(26) 



(27) 



in any consistent theory py\. Note that the Steinberger calculation is then a 
special case wherein N c = 1 and this feature is the reason that his calculation 
agrees with the usual quark model result. (Of course, if one also examines 
other anomalous processes such as f]° — > 7r + 7r~7 agreement between experi- 
ment and theory does require N c = 3, and this provides the real proof.) 

We can now check the validity of the various conservation conditions — Eq. 
[T5l the first of which reads 



QiT„ vl {qi,q2) 



+ 



le 

T 
i 



d 4 S 

,1 



Tr 



-)7 7 7s 



(28) 



However, the integrals which involve a single factor of photon momentum qi 
or §2 vanish, since the epsilon tensor associated with the trace 

^l^lulalfilb = 4ze A1 „ Q/3 , 

requires contraction with two independent momenta in order to be nonzero. 
Thus, defining 

W^(s) = Tr (Ly v , ) ^ 7 7 7 5 ) (29) 
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we have 

ic^ 1 r cfc^s 

?fWgi,? 2 ) = ~2j j^ w ^ s + ?0 " + 92)] (30) 

If the integrals in Eq. [30] were convergent, or diverged no worse than loga- 
rithmically, then we could shift the integration variables freely and thereby 
obtaining zero and verifying gauge invariance. However, because there ex- 
ists a linear divergence at large s we must be more careful. Using Taylor's 
theorem 

/ ^Y F{s + a) = I^Y [F{s) + aadaF{s) + • • - ] (31) 



we can write 



<£W<&,fc) = -^fe - / -^y\d a W ul {s) + ...} (32) 

where the higher order terms, denoted by the ellipsis, vanish, while the piece 
we have retained can be evaluated via Gauss' theorem, yielding 



Similarly we find 



and 



qiT^(qi,q 2 ) = -^-^ysqi q 2 (33) 



<l2T lun (q 1 ,q 2 ) = —e lunS <fi<fi (34) 



(9i + Q2) J T l , ul (q 1 , q 2 ) = (35) 

Thus gauge invariance is violated, which would have serious consequences for 
photon interactions. This problem can be solved by appending a polynomial 
in the external momenta, which can be done without affecting the absorptive 
component of the amplitude. Thus defining the physical decay amplitude via 

we see that gauge invariance is restored — 

q^ s (qi, fc) = q 2 t^{q u q 2 ) = (37) 
21 



However, the axial contraction now yields 



ie 2 



fa + Q2)X h A s (^ <&) = 4^3 (38) 

so that the axial current is no longer conserved. Thus, the axial symmetry has 
been broken via proper quantization of the theory — there exists an anomaly. 
We have then 

C"(«)=^^ (39) 
which corresponds to the operator condition 

2 



Using the PCAC condition, we have then 
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< 2 7 |d"jf M |0 >= F v ml— < 2 7 |vr° >= e^q^q^ (41) 

Of course, the above evaluation is a simple LO perturbative calculation 
and it is interesting to ask what is the effect of interactions, whether these 
be higher order electromagnetic or strong interaction effects. The answer, 
as shown by Adler and Bardeen[Bn], is that such effects do not modify the 
chiral anomaly. The basic argument is that any such interactions will result in 
changes to the triangle diagrams in which the linear divergence is removed — 
the diagrams become more convergent. Because of this modification the 
vanishing of Eq. [SU] from such diagrams is assured and we conclude that the 
lowest order calculation given above must be true to all orders. 



3.3 Alternative Approaches: All Roads Lead to Rome 

What are some of the alternative methods that can be used in order to deal 
with the short distance behavior? Some of the possibilities include: 

i) Pauli-Villars: As has already been mentioned, a Pauli-Villars regulator 
can be used in order to make the results finite. In this procedure 
one defines the physical amplitude as the difference of the amplitude 
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calculated as above and the same amplitude calculated with quarks 
with large mass M — 



T phymcai^ ^ = ^ T^xfai, ft) - T^ x (q u ft) (42) 

where T^ A (gi,ft) is identical to T^^i, g 2 ) but with the massless 
fermion propagators replaced by propagators having mass M. That 
is, 

Z£y(«i, 92) = ^ 7 (ft, ft) + 0<?2, 9i) ( 43 ) 

where 

ie 2 r d 4 s „ / 1 1 1 



^ 7 (ft,ft) = ^,Tr | v /+ ^_ M 7, 7 f ]? 7 v _^_ M 7 7 7 5 

(44) 



Taking the axial divergence we find 



-,2 



<h) X U£ x (<li, ft) = Wflf 9? (45) 



e 2 



so that 

(gi + g 3 ) A TjS* lca, (?i > ft) = ^Kpe^q^q^ (46) 
We reproduce in this way then the axial anomaly 

16tt 5 



5V 5 3 A = (47) 



ii) Path Integration: Path integral methods can also be used. In this case 
one quantizes by using the generating functional, which sums over all 
field configurations 

W = J [#] [#] [dAp] exp iSty, Ap] (48) 

Now under a local field transformation such as Eq. [12] or Eq. [13] the 
action becomes — 

Sty', ft, AJ = Sty, A,} + / d 4 x(3(x)d» JJx) (49) 



Fujikawa showed that under a vector transformation — Eq. [12] — the 
integration measure is unchanged — 

[#'] [#'] [<*4J = [#] [#] [d^] (50) 



23 



so that, requiring equality of the two representations for arbitrary f3(x) 
yields 

= (51) 

i.e., the classical field symmetry is also a quantum field symmetry. 
However, this is not not the case for a local axial transformation — Eq. 
IT51 In this case we have a non-unit Jacobian 

Wa\ Wa\ [dAJ = [#] [#] [dA„] J (52) 

where 

J = exp(— 2iTr7(x)75) (53) 

Here the trace is over not only the Dirac indices but also over spacetime 
and must be regulated in order not to diverge. Employing a covariant 
regulator of the form 

exp(- l/) 2 /M 2 ) 

which cuts off the high energy (short distance) modes, and the result 



V> = D 2 + -j-^F^ (54) 
Fujikawa demonstrated that 

2 

J = expH / d 4 x 7 (x)-^F^i>0 (55) 

The connection with the chiral anomaly can be made by noting that 
the generating functional 



W = J [dip] [#] [dAp] exp iS[if>, V, A„] (56) 
after a local axial transformation Eq. [131 assumes the form 



e 



2 



SW,i>,A„] + / d A X1 {x){d x JUx) - —F, u Fn 

(57) 

so that invariance of the functional integration for arbitrary 7(2;) yields 
the anomaly condition 

d X JL = 4^ F ^ U (58) 
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Point splitting: Since the problems arise when the field and its conju- 
gate are at the same spacetime point, one can define the axial current 
via the definition 



1 , 1 



1 



^(x) =lim-ip{x + -e)-T 3 'y^ 5 'ifj(x--e)exp{ie J^^ dy p A p {y)) (59) 
If we now take the divergence we find 

k 1 - 1 



lim 



[X 



2 e) 2 T315 



e4(x + -e) 



e# - if) 



ij){x - -e) 



1 1 1 l r x+ \ e 

ei/j(x + -e)-r 3 ^ 5 ip(x - -e)e v d*A u j exp(ie J dy p A p {y)) 

\ - 1 1 1 /• a; +5 e fl 

limee^ ^{x + -e)~T 3 7 M 7 5 ^(x - -e)F^exp(ie / i dy p A p (y)) 



In taking the limit e — > we use the short distance behavior of the 
Dirac field — 



limTr 

e->0 



1 1 1 > 

e^-TaT^Gtx - -e,x + -e, 



167T 



1 P v 



to once again we obtain the axial anomaly 

„2 



«9 A J, 3 



A 



T?^ v TP 

167T 2 



(61) 



(62) 



(60) 



) Geometric Approach: Using geometric methods, Bardeen was able to 
identify the full form of the anomaly|61j. In the case of SU(2) and 
electromagnetism the anomalous Lagrangian density is 



N K 



48vr 2 
with, defining 



eA^Ti{QL u L a Lp - QR v R a Rp) + ie 2 F lw A a T^ 



(63) 



U = exp(fr • (^jt/Ft,. 



duUU\ 



Tp = Ti (Q 2 Lp-Q 2 R p + ^QUQtfLp-^QU^QURp) 
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The piece of Eq. [2H responsible for 7r° — > 77 can be found by expanding 
to first order in the pion field 

p 2 N 

£a = T | z ^Tr(QV 3 )e^F^ Q <V° 



F^n (65) 



e 2 N ( 



c 



2An 2 F 7r 

which once again reproduces the anomaly prediction. 

In each case then one is forced to modify short distance properties in 
order to produce a consistent quantum field theory, and it is this feature 
which breaks the classical symmetry and produces the anomaly. No matter 
how it is obtained, the result in the case of the two photon decay amplitude 
of the neutral pion is that the decay amplitude is precisely predicted to be 

e 2 a 

^7r°77 = ^2p e ^al3 e l e 2 ( ll ( l2 (66) 

leading to a decay rate 

2 S 
OL 771 

r ^ = 64^ = 7 - ?6eV (6?) 

7T 

where a = e 2 /4ir is the fine structure constant, in excellent agreement with 
the experimental value. However, before a careful confrontation with exper- 
iment, it is essential to consider corrections which arise from the the feature 
that the prediction Eq. [67] is made in the limit of chiral symmnetry, rather 
than the real world. 



3.4 Real World Corrections 

As mentioned above, the prediction Eq. [H7J is unsatisfactory in that the 
chiral limit, in which both the u, d quarks and the pion are massless does 
not represent the real world. The u, d quarks are light but not massless, and 
in turn the pion is the lightest hadron but certainly does not possess zero 
mass. More importantly, because the light quarks are nondegenerate, the 
physical tt° meson is not a pure U(3) state \P% > but is instead a mixture 
of \Ps >, \P$ >, and \Pq > states. This mixing is important since a simple 
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U(3) picture of the decay predicts the bare amplitudes for the two photon 
decay to be 




Ap 311 : A Pg77 : A Poll -l:\-:2\- (68) 



so that, even if the mixture of \P$ >, \Pq > states in the neutral pion is 
small, they can play an important role in the ir° — > 77 decay amplitude. It is 
important to note here that the predictions given in Eq. [68] are given U(3), 
wherein the wavefunctions of the 7]$ and 770 are taken to be identical. This 
is perhaps surprising, since the rj 8 is a Goldstone boson and is massless in 
the chiral symmetric limit, while rj is not. Its mass is known to arise due to 
the axial anomaly, but does vanish in the N c — > 00 limit [13]. Nevertheless, 
here and in other circumstances the use of U(3) symmetry is known to give 
good results [62J. In terms of chiral perturbation theory, the prediction of the 
chiral anomaly for the tt° — > 77 decay amplitude, since it involves the four- 
dimensional Levi-Civita tensor, is already four-derivative-(9(g 4 ) — and is two 
orders higher than the leading order strong interaction chiral effects, which 
are C(g 2 )[63]. Modifications of the pion decay amplitude due to particle 
mixing are also 0{q i ). Since the pion is isovector and the 77 and 77 are 
isoscalar the mixing matrix elements must be proportional to 

m u - m d . 

Of course, there are in addition corrections of 0(q 6 ) which involve factors 
such as m^/A 2 , but these are presumably higher order and smaller than those 
which come from mixing. 

It is useful to make a simple back of the envelope calculation of the 
modifications due to mixing effects by use of the representation of the physical 
states |7r°,?7,?7 > in terms of the bare U(3) states P = |7r 3 , 7/8,770 >. In order 
to do this it is useful to first review the P — > 77 decay rates predicted by 
the anomaly 

A Prr = T^FT (69) 



16tt 2 F p 

r 



-^^77 I p 



7T77 



Itt 



where m P ,F P represents the mass and decay constants of the pseudoscalars. 
In Eq. [68] the relative amplitudes of the amplitudes was given in lowest 
order, where all of the decay constants are equal. In the next chiral order 
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there are corrections to thisjOI] which lead to F 8 /F 3 = F v /F n ~ 1.25 and 
F /F 3 = F 1 1 F n ~ 1. We shall use these values and amend the values 
in Eqj68] by them. With this the two gamma widths are calculated to be 
T(r) — > 77) = 0.11 keV and T(r] — > 77) = 7.40 keV. These are not in 
agreement with the experimental results [H] T(r) — > 77) ~ (0.51 ± 0.05)keV 
and T(?7 — > 77) ~ 4.28 ± 0.38 keV. The fact that the 77 calculated width is 
too low, and the 77 too high, is strong evidence for 77, 77 mixing. For this one 
can writel 



\r] > = cos^^g > — sin^o > 
1 77 > = sin^^g > + cos^^o > 



(70) 



Using this representation one obtains for the decay amplitudes 



.4 



a 



A/ 77 



a 

4vrV3 



cos 9 sin^v^ 



F 1 



(71) 



sin 6 cos^y^ 



F 1 



Treating 9 as a free parameter one can obtain agreement with the experi- 
mental values [T4~j with 9 between -20 and -25 degrees [64"]. 

There is another method to calculate the 77, 77 mixing which is instructive 
and will be helpful in estimating the chiral corrections to the 7r° decay rate. 
For the 77, 77 mixing amplitudes one diagonalizes the mass eigenvalue matrix 
to obtain the physical 77, 77 states 



M 2 



m 



ON 



m 



OS 



in?. 



M 2 








where m 8 ,m m,m s, the masses of the unperturbed 778,770 states and the 
off-diagonal mass mixing matrix element which we can take as a parameter. 
The eigenvalue equation is in terms of the square of the masses since it is for 
Bosons. There are two eigenvalues for M which are m v ,m >. To solve this 
eigenvalue equation we need the value of m 8 which can be obtained from the 
relationship 



{Am 2 K ~ ml){\ + 5) 



1 

3 

1 2 1 2 
-{m K0 + m K , 



(72) 
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For 6 = this is the Gell-Mann-Okubu mass formula which is derived on the 
basis of SU(3). Since this is not exact there are chiral corrections for which 
the leading log estimate gives 5 ~ 0.16164] . The value of mo can be obtained 



by observing that from the trace of the mass matrix in Eq l3.4l m§ + m\ = 
m 2 + m 2 / . Taking 5 as a free parameter the eigenvectors (specified by the 

angle 9 of Eq JTT]) can be obtained and from this the values of T(r), rj' — >)7, 7). 
The data can be fit in the range of 5 from ~ 0.16 to 0.22 in approximate 
agreement with the leading log estimate [64]. 

Now we are ready to estimate the magnitude of rj, rj mixing in the ir° 
amplitude, which can be approximately written as 

|tt° >~ |tt 3 > +9 ri \rj > +9 V > \rj > (73) 

The mixing amplitudes 0^,9 ' can be obtained using perturbation theory. 
The off diagonal matrix elements of the mass matrix(see e.g. [13]) and the 
mixing amplitudes in terms of the quark masses are 

m 38 = < V8\M\tt 3 >= -J=(m d - m u ) 
m 30 = < Vo\M\ir 3 >= - m u ) 



9 



cos #m| 8 — sin 6*m| 

v 9 2 

m^ - m^o 



q , = sinflmls + cos 9m 2 30 
T ' m 2 , — m 2 

r? 7r u 

where B is related to the vacuum expectation value of the quark scalar 
density via 

m 2 

£0 = ~ T , < 0|w|0 >^ ^ (75) 

Using the value of Bq in terms of m n o and using the quark mass ratio 
m u/ m d — 0.56 [T7] numerical calculations can be performed. A number of 
interesting features emerge. First the decay rate is increased by ~ 4%. This 
value is only weakly dependent on the parameter 6 which corrects the GMO 
mass formula EqJ73l For 5 ~ 0.18, for which the 77, rj rates are in agreement 
with experiment, the values of 9 V ~ 0.015 rad and 9 > ~ 0.0032 rad. The 
value of r(7r° — > 77) ~8.1eV a ~ 4.5% increase over the value predicted by 
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the chiral anomaly. The contribution of the 77 is ~ 3 % and from the 77' ~ 
1%. As we shall see, it is remarkable that this back of the envelope estimate 
is in agreement with chiral perturbation theory calculations. 

Stimulated by the PrimEx experiment, calculation of QCD corrections 
to the chiral anomaly prediction for the ir° — > 77 decay amplitude has been 
undertaken by a number of groups with remarkably similar results: 

a) An estimate by Ioffe and Oganesian|65] using sum rules but including 
only 7T°, r]° mixing yielded a 3% enhancement — r„-o^ 77 = 7.93 ± 0.12 
eV. However, in view of the mixing discussion given above, the neglect 
of 7/ effects is obviously a serious concern. 

Other authors have included all corrections within the context of various 
approximations to QCD: 

b) The work of Goity, Bernstein, and Holstein[29] involves the use of chiral 
£7(3) x U(3) symmetry and the large N c limit in order to include the 
77' as a Goldstone boson and includes all corrections to the anomaly 
prediction of 0(q 6 ) and 0(q 4 x -±). These corrections predict a 4.5% 
enhancement to the decay rate — F^o^.^ = 8.13 ± 0.08 eV. 

c) Another approach was that by Ananthanarayan and Moussallam|31] 
who looked at chiral perturbation theory in the anomaly sector with 
the inclusion of dynamical photons. In this way they looked both at 
quark mass effects and at electromagnetic corrections of order q 4 e 2 . 
The result was a predicted decay rate — r^o^^ = 8.06 ±0.02 ±0.06 eV, 
in excellent agreement with the Goity et al. calculation. 

d) An alternative approach by Kampf and Moussallam[30] involved the 
use of chiral perturbation theory but with a modified counting scheme 
in which m s is considered of order q. The result was a prediction — 
rv^ 77 = 8.09 ±0.11 eV. 

These values are plotted in Figj5j There is obviously remarkably little 
scatter among these theoretical calculations that chiral symmetry breaking 
quark mass effects should increase the decay rate of the neutral pion from 
its 7.76 eV value predicted by the anomaly to 8.10 eV, the average of these 
results. The basic reason responsible for this 4.5% enhancement can be 
found in the pseudoscalar mixing estimate given above. However, in order 
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to confirm this enhancement it is clearly necessary to perform an experiment 
looking at ir° — > 77 decay at the percent level. On the theoretical side, 
the possible ~ 1% errors in the 7r°77 decay rate estimates arise not from 
convergence issues in the chiral expansion but rather from uncertainty in 
some of the calculational input parameters as well as the mixing estimates 
due to isospin breaking. However, since these predictions are already at 
the one loop level in the chiral expansion, it is unlikely that they will be 
significantly improved in the near future. 

3.5 Other Anomalous Processes 

From the form of the full SU(2) anomalous effective Lagrangian Eq. (63] it 
is clear that, in addition to the triangle diagrams which we have considered 
above, there also exist anomalous box and pentagon diagrams. The presence 
of the Levi-Civita tensor means that the processes described by Eq. [63] 
are of a different character than those described by the conventional chiral 
Lagrangian (63]. Witten has pointed out that under what he calls an "intrinsic 
parity transformation" under which the pseudoscalar fields undergo a change 
of sign but spacetime coordinates remain unchanged |66j- 

8 ' 8 

<t> P ^-<t> P so f/ = exp(^£A^f)^exp(^£A^f) = t/t 

tp 3=1 tp 3=1 

— the conventional chiral Lagrangian ala Gasser and Leutwyler |63] . [67] . which 
contains terms such as 

Trd^d^U, [Tro^f/tfW] 2 , Tr(d fl U^d u U)^(d^d u U), etc., 

includes only the interactions of an even number of pseudoscalars (or axial 
currents) or the coupling of one or two photons to two, four, six, etc. pseu- 
doscalars or axial currents, while the SU(2) anomalous Lagrangian of Eq. [63] 
(and its SU(3) generalization) is symmetric under U —> U' and L M — > 
and so describes transitions between processes involving an odd number of 
pseudoscalars (or axial currents) or of the coupling of one or two photons 
to one, three, five, etc. pseudoscalars or axial currents |66j. The existence of 
the chiral anomaly then, in addition to predicting the decay amplitude for 
n — > 77, also makes no-free-parameter predictions for processes such as 

a) 77T° — > 7T + 7T~ 
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r(rc -> yy) Chiral Calculations 
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Figure 5: The r(7r° — > 77) width in eV predicted by the chiral 
calculations |29 [ |30 | I3T] . The lower dashed line is the predictions of the chiral 
anomaly[3j g] (r(7r° ->■ 77) = 7.760eV, t(tt ) = 0.838 • 10~ 16 sec). The upper 
solid line is the chiral prediction and the dotted lines show the estimated 1% 
error[2HJ EUJ EI](r(7r -> 77) = 8.10eV, r(7r°) = 0.80 • 10~ 16 sec). 
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b) 7r + — > e + u e '~f 

c) etc. 

or, by extending our analysis to the case of SU(3), for reactions such as 

d) K+K~ -> tt+tt-tt 

e) 1] — > 7r + 7r _ 7 

f) K + -)■ n + n~e + v e 

g) K + -> e + z/ e 7 

h) etc. 

While it is possible in principle to probe the validity of the chiral anomaly 
by measurements of any these reactions, analysis of the decay (f) gen- 
erally assumes the value for the appropriate vector form factor predicted by 
the anomaly. In the radiative 77 decay reaction (e), significant mixing with 
the 7/ obscures a precise test of the chiral anomaly |68j. In the case of the 
737T vertex (a) a test was performed many years ago by Antipov et al. with 
the result ESI, E3 



Amp^ = (12.9 ± 0.9 ± 0.5) GeV 3 vs. Amp—" = TT^k = 9 - 7 Gey3 

12.71 r n 

(76) 

Obviously, an up to date modern measurement of the 737r amplitude would 
be of interest. 

A second reaction that has been utilized as a test of the chiral anomaly is 
that of radiative pion decay (b). In fact, the pi-beta experiment mentioned in 
the introduction uses a measured weak polar vector form factor to predict the 
neutral pion decay amplitude by isospin invariance and is used to generate 
the PDG average[2"5] . Because of this fact, we discuss here this process in a 
bit more detail. The transition amplitude can be written in general as 

M^e+v^ = -^=-V ud M^(p, qy*u(p u )Y{l + l5)v{p e ) (77) 
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where Gp is the weak decay constant, V u d is the CKW parameter connecting 
the u and d quarks, and[T3"] 



M^PA) = J <Fxe** x <0\T(j; m (x)J^ 2 (0)\n + ((p)> 

= -V2^ (p iy_ m 2 ) < 7r+(p - q\J^ + {p> +V2F^ 

- h A ((p - q)^q u - rj^q ■ (p - q)) + ih v e iWOL pq OL pP (78) 

Here the second line represents the Born diagram together with a contact 
term required for gauge invariance and in the following line the subscripts V 
and A indicate whether the weak vector or axial-vector of the weak currents 
are involved. The purpose of the pi-beta experiment was to measure the rate 
of the pion beta decay reaction ir + — > 7r°e + b> e as a test of the conserved vector 
current (CVC) relation or, by using CVC, to measure the CKM parameter V u d 
in a reaction where strong interaction uncertainties are not relevant. Because 
this is such a small (~ 10 -8 ) branching ratio process, the experiment also 
detected other higher branching ratio reactions such as radiative pion decay 
7r + — 7- e + + v e + 7. Ordinarily a radiative decay is primarily sensitive to the 
Born amplitude, which simply generates a correction to the non-radiative 
decay process via 

dT a dT Q 
~ + . . . 

dVtdk k d£l 

so that direct decay amplitudes are hidden under this huge bremsstrahlung 
background. However, because the nonradiative decay process tt + — > e + u e 
is helicity-suppressed the direct decay amplitudes hy, h A above can both 
be measured and this was done by the pi-beta experimenters. The direct 
axial-vector decay amplitude h A determines, via PCAC, the charged pion 
polarizability, while the direct polar-vector amplitude hy is related to the 
pion decay constant via|13] 



2 e 2 /1 'r°77 



That this should be the case is clear from Bardeen anomaly Lagrangian Eq. 
(k?1 or simply from isotopic spin invariance. Indeed a simple isotopic spin 
rotation relates the the neutral pion decay amplitude 

AW^i) = e 2 / d 4 xe^ x < 0\T(j; m (x)Jr(0))k°(p) > (79) 
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to the direct polar vector amplitude hy defined in Eq. [78] via a Clebsch- 
Gordan coefficient. Thus a measurement of hy can be used to give an ex- 
perimental value for the ir° — > 77 decay rate. The result of the pi-beta 
experiment yields [25] 

rpr^ eta = 7.7 ± 1.0 eV, 

and this number is used as an input to the PDG average. (Note, however, 
that since isospin invariance is broken at the 1% level, the present 13% 
precision of this method is not an issue, but the uncertainty associated with 
isospin violation ultimately limits its use at the 1% level unless this breaking 
is understood theoretically.) 

While discussion of any of the processes listed above and their connection 
to the chiral anomaly is certainly of interest, in the present paper, we shall 
focus on the ir° — > 77 process. 

3.6 Physics of the Anomaly 

Above we have shown how the chiral anomaly leads to a remarkable successful 
agreement between experiment and theory for the decay rate of the neutral 
pion. However, this derivation is somewhat formal and it remains to show 
what the "physics" of the anomaly is — that is, why must quantization destroy 
the classical axial symmetry. In order to present an answer to this question 
it is useful to first examine the Schwinger model, which is the name usual 
given to massless electrodynamics in one plus one dimensions |54j. Here the 
Lagrange density is given by 

C = ^ fl^ - -F^F» V (80) 

where 

D li = d ll + ieA lt (81) 

is the covariant derivative and the 2x2 Dirac matrices are given in terms 
of the Pauli matrices via 

7 = 0-1 and 7 1 = ia 2 (82) 
It is then easy to see at the classical level that we have equations of motion 

i pip = and UA^ = ej M (83) 



35 



where 

3n = (84) 

is the vector current and is conserved — <9 M j At = 0. It is easy to see that there 
also exists an axial current 

jl = ^lifl^ with 75 = -7°7 1 = o- 3 (85) 

which is conserved — d^j^ = 0. For later use, we note also that 

3l = (86) 

where e^ u is the two-dimensional Levi-Civita tensor. So far then, this looks 
simply like a two-dimensional version of massless QED. 

However, upon quantization it can be shown that the Lagrangian can be 
written as 

C = $i ft, - \f, u F^ - (87) 

that is, in terms of a noninteracting system of massless spin 1/2 particles 
and free "photons" having mass m 2 = e 2 /ir. Also in the quantized theory 
the axial current is no longer conserved. Rather we have 

Pjl = -%f F *» ( 88 ) 

That is, axial current conservation is broken by quantization — there exists 
an anomaly. 

The physical origin of the anomaly can be seen via an argument due to 
Widom and Srivastava[7T] by considering the vacuum state in the quantized 
theory, which, according to Dirac, can be considered as a filled set of negative 
energy states. In the absence of an external electric field there exists a density 
of electron states dp/2-ir with momenta evenly distributed between p = —00 
and p = 00, and so there is no net current. Now consider what happens 
in the presence of a constant electric field E — a net current flow develops, 
which increases with time. In terms of the current density j we have 

dj r°° dp dv 

dt ~ e J-oo 2vr ~dl ^ ' 

where, using the Lorentz force law — dp/dt = eE — we have 

dv d d 1 — eEm 2 . . 

— = — — Jm 2 + p 2 = r 90 

dt dtdp y ( m 2 +p 2^ 



1 2 



3(3 



Performing the integration we find a result 

dj e 2 E . . 

—r = 91 

dt it ( J 

which is independent of the mass. Since the vacuum charge density A is 

independent of position — dX/dx = — we have, defining j M = (A, j) and 
using Eq. [86] 
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which is the chiral anomaly. Note that Eq. [92] can also be written as 

= (ef + ^A"J (93) 
In Lorentz gauge — d^j^ = — we have 

U = ( 94 ) 
so that the equation of motion in Eq. [831 becomes 

(□ + -)A, = (95) 

which indicates that the "photon" has developed a mass m 2 = e 2 /ir. We see 
then that in this picture the origin of the anomaly is clear and arises from 
the feature that the vacuum state of the quantized system is altered in the 
presence of an applied electric field. 

An alternative way to understand the same result has been given by 
Jackiw[72"] . wherein one looks at solutions of the time-independent Dirac 
equation 

Eip = 7 o7i f-i^j; - eAj (96) 

In the case of a constant vector potential then there are two classes of solu- 
tions 

i/}+(x) = I j with E = p — eA 

•*!>-{*) = ( JL) with E = -p + eA (97) 



37 



E 




Figure 6: The Dirac Sea. 

where the subscript ± specifies the chirality of the solution — i.e., the eigen- 
values of the operators |(1 ± 75). If A = we see then that the vacuum 
consists of (negative energy) states with p < for positive chirality and 
p > for negative chirality — cf. Figure [6j 

Now suppose we make an adiabatic change from A = to a nonzero field 
with A = e. In the presence of the field, the vacuum states become those 
with p < ee for positive chirality and p > ee for negative chirality, meaning 
that there is a net chirality production 

r ee dp ee , s 

A X = 2/ £ = - (98) 

JO Z7T 7T 

This result should be expected from the chiral anomaly, which requires a 
time-varying axial charge 

d „ e ,_, e dA . . 

-Q 5 = -E=-— 99 
at n n dt 

Since 

Q 5 = J (fa^W (100) 

corresponds to chiral charge, we find, integrating both sides of Eq. 1991 

e pf 

A X = —AA = — (101) 

7T 7T 

in agreement with Eq. [98j Once again we see that the origin of the anomaly 
is the modification of the vacuum in the presence of an applied electric field. 
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A simple handwaving argument can be used to generalize the latter ar- 
gument to four spacetime dimensions [73]. A constant magnetic field in the 
^-direction can be represented by the vector potential 

A=-fxB (102) 
The energies in the presence of the field are given by the Landau levels 



E n = ±yJp 2 z + 2eBk with n = 0,1,2,... (103) 

If we now turn on an electric field, also in the z-direction, then the energy 
levels change in accord with the Lorentz force law 

^ = eE (104) 
dt v ; 

As p z changes the energy levels change, but negative energy levels remain 
negative so that no particle creation occurs. However, if k — then the 
structure of the levels is as in the case of the Schwinger model, so that 
helicity is produced at the rate 

where p is magnetic flux density. Since magnetic flux density is quantized in 
terms of 27r/e, we have 

„=§ (106) 

and thus 

Since helicity and chirality are the same for a massless system we can write 
this as the covariant equation 

which is the standard chiral anomaly. 
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3.7 Theoretical Summary 



We have seen above that a consistent quantum field theory requires modifi- 
cation of the naive short distance behavior and leads to the breaking of axial 
symmetry — 

This anomalous symmetry breaking leads to a decay rate 

T™ = 7.76 eV 

When real world corrections such as chiral symmetry breaking and mixing 
are included this prediction is raised by about 4% to 

rf4° 77 = 8.10 eV, 

with remarkably little theoretical uncertainty, and a (9(1%) experimental 
verification of this prediction is also a validation of QCD. 

Although rigorous quantum field theoretical arguments lead unambigu- 
ously to this violation of axial symmetry it is also useful to understand the 
"physics" of this phenomenon and we presented arguments which show that 
origin of the anomaly is the modification of the vacuum state of the field 
theory in the presence of an electromagnetic field. 

Having then understood the connection of the decay rate of the ir° with 
QCD, we move now to examine the experiments, on which the PDB number 
is based. 



4 Measurement of the tt Lifetime: The Par- 
ticle Data Book 

As discussed in section [21 by the end of 1965 it was known that r(7r°) ~ 10- 16 
sec, based on measurement of the 7r° decay distance[2T] and determination 
of the Primakoff cross sectionjIT]. This section will cover experiments per- 
formed during the period from 1970 through 1988, on which the present 
value in the particle data book is primarily based [14]. As an overview we 
presented in FigJT] the results which will be presented in this section, the 
recent Primakoff measurement [27J along with the current theoretical predic- 
tions. The first three measurements were performed using the Primakoff 
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effect [2SJ E3J E2] at Frascati, DESY, Tomsk, and Cornell respectively. The 
fourth result is a direct measurement of the 7r° decay distance [23] performed 
at CERN at high energies, and the fifth result was obtained via a ir° produc- 
tion cross section measurement in e + e~ collisions performed at DESY [21]. 
The sixth result is due to a measurement was performed of radiative pion 
decay(7r + — > e + i"y)^I5\. Using isospin invariance, the weak polar- vector form 
factor contributing to this decay channel is related by a simple isospin rota- 
tion to the amplitude for 7r° — > 77(see Sec J3.5l for further discussion). The 
last point is a recent Primakoff effect measurement |27] (discussed in Sec. 5). 
We now discuss each measurement that is used in the present version of the 
Particle Data Book[HJ [2Sj and present an overall outlook. 

4.1 Primakoff Effect Measurements 

As mentioned above, in the decade between 1965 and 1974 there were four 
experiments performed [4T] |2"5] 153] E2] which used the Primakoff effect [22J. 
All of these early experiments utilized bremsstrahlung photons, which have 
continuous energies below E , the energy of the electrons which produced 
them. A measure of the photon energy was obtained from the opening angle 
distribution from the two photon decay(see the appendix at the end of this 
section). The Primakoff experiments determine the cross sections for the 
7 + A — > 7T° + A reaction. At small angles this reaction is dominated by the 
7 + 7* — > 7T° process, in which one of the gamma rays is due to the Coulomb 
field of the nucleus, which remains in its ground state. However, this means 
that the neutral pions are also produced by the photons interacting with the 
nucleons and either leave the nucleus in its ground state (nuclear coherent 
production) or leave the nucleus in excited or continuum states with the 
production of pions (nuclear incoherent). The Primakoff effect dominates 
at very small angles, but the small contributions of the nuclear effects must 
be subtracted from this small angle signal. In practice this is accomplished 
by fitting the parameters of model calculations for the nuclear effects using 
the larger angle data. As the accuracy of the experiments improved, these 
requirements have correspondingly become more stringent. We shall now 
sketch the basics of the Primakoff and nuclear effects before examining the 
experiments. 

The photoproduction of pions from a complex nucleus, 7 + A — > ir° + A, 
can be described by the sum of Coulomb T c and Strong T s amplitudes. 
Including incoherent production, the differential cross section is: 
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T c | = T 77 — ——\F e . m \Q)\ sm ff 




Q 2 



|r s | 2 =c 5( T^ 2 |F s (g)| 2 sin 2 ^ 

(m 2 /2E) 2 + E 2 sin 2 W 



(109) 



where Tc, Ts are the Coulomb (Primakoff ) and strong amplitudes. The phase 
cf) originates from the 'jp — > ir°p amplitude and is fitted to the data. The first 
two terms in the first line represent the coherent cross section for which the 
nucleus is left in its ground state. dai/dn{i = P, S, inter, inc) are the cross 
sections for Primakoff, strong, interference, and incoherent processes (the 
latter involving target nucleus excitation or break up). T 77 is the pion decay 
width (the primary objective of the Primakoff experiments), Z is the atomic 
number, m, (3, 6^ are the mass, velocity and production angle of the pion, 
and E is the energy of incoming photon. For the spin targets employed 
in these experiments the coherent cross sections are non-spin flip and as a 
consequence have a sin 2 6^ dependence, ensuring that no scattering occurs 
at forward or backward angles. Q is the momentum transfer to the nucleus, 
and E e m (Q), Fg(Q) are the nuclear electromagnetic and strong form factors, 
corrected for final state interactions of the outgoing pion [711 [751, ITo] . (Note 
that due to the absorption of the outgoing pions these form factors are com- 
plex.) The shape of the strong cross section das/dn is determined by the 
dependence of the absolute value of the strong form factor | F$(Q) | and the 
sm(6 n ) 2 factor, a^ is the the non-spin flip part of the neutral pion photopro- 
duction cross section on the nucleon which the Cornell group parameterized 
as <t at ~ l00E 2 fib^52\. The experimental results are fit with the theoretical 
cross sections with four free parameters T 77 , Cs, Ci nc , <fi. The fitting parame- 
ter Ci nc , which is not shown in Eq. 11091 is introduced to vary the magnitude 
of the theoretical incoherent cross sections. This latter contribution is small. 

The Primakoff cross section peaks at small angles and is quite narrow. 
The general features of the cross section can be seen by following the treat- 
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ment of Gourdin[77] in the high energy limit which reads 



jn ^77 m 3 f\2 \Fe.m.(Q)\ f(Q /Q m in) 

fit = Q 2 /Q 2 mm ) = (l-l/t)/t 
QLn - (rn 2 /2E) 2 9 P:max ~ m 2 /(2E 2 ) (110) 

where m is the mass of the produced meson. The energy-independent func- 
tion f(t) is shown in Fig. 11101 and can be seen to rise rapidly from forward 
angles (t= 1) and peak at t = 2. From this behavior we can easily infer the 
angle at which the Primakoff cross section is a maximum — 6p :max shown in 
Fig. The numerical values are plotted as a function of photon energy. 
The small value of 6p- max means that the experiments that can measure the 
Primakoff effect must detect 7r° production at small angles. In addition, from 
the shape of f(t) it can be inferred that the width of the peak is ~ 29p- max so 
that the detector needs to have excellent angular resolution. Furthermore, in 
order to suppress background, it is very useful for the detector to have good 
energy resolution as well. 

The general features of the cross section for 7r° production are shown 
in Fig. |HJ Here the contribution of the Primakoff, nuclear coherent, and 
their interference cross sections are shown (the small contribution of the 
nuclear incoherent cross section is not given). The figure demonstrates the 
dramatic increase in the cross section as the photon energy is increased and 
also indicates how the Primakoff cross section increases relative to the nuclear 
coherent cross section for a heavy nucleus relative to a lighter one. The 
reason for this increase is nuclear absorption of the outgoing neutral pions. 
If this effect were absent the nuclear coherent cross section would scale as 
A 2 , where A is the nuclear mass number. However, for strong absorption 
the cross section only increases as A 2 / 3 , which is close to the actual physical 
case. 

The early Primakoff experiments were performed at Frascati in 1965 [4"T] 
using bremsstrahlung photon beams with end point energy E ~ 1 GeV and 
a Pb target. The experiment was then repeated at DESY in 1970 by a similar 
group with Eq = 1.5 on 2.0 GeV and C,Zn, Ag, and Pb targets [28]. Another 
experiment was performed at Tomsk in 1970 with Eq = 1.1 GeV and a Pb 
target [33]. These experiments were analyzed with the early nuclear coherent 
calculations of Morpurgo|74j. This calculation was a significant advance over 
previous ones which neglected final state interactions of the outgoing pions 
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and explained the main features of the nuclear coherent photoproduction, 
which included absorption of the outgoing pions in the nuclear medium. 
However, the further approximation of a uniform nuclear density was made, 
which limited the accuracy of these calculations. In addition it was later 
shown that they did not include the effect of rescattering of the outgoing 
pions back to small angles [75]. In addition, the phase of the interference 
amplitude was assumed to be angle-independent. This phase is important 
since this amplitude is the only correction that must be applied to the small 
angle cross section in an accurate determination of the r(7r° — > 77) width. 
We therefore conclude that, due to the low energy of the first three Primakoff 
experiments and the use of the Morpurgo calculation [71], these experiments 
are not sufficiently precise to include in a modern average of the ir° lifetime. 
This conclusion differs from the current particle data book average [14] which 
includes both the CERN[28] and Tomsk [33] measurements. 

The last of the early Primakoff experiments was performed at Cornell 
using bremsstrahlung beams of endpoint energy E = 4.4 and 6.6 GeV with 
Be,Al,Cu,Ag, and U targets [52]. This experiment had the advantages of 
higher beam energies and also utilized the improved calculations of Faldt [75] . 
This calculation was the first to include the effect of rescattering of the out- 
going pions back to small angles and also included the real part of pion 
rescattering in the final state. These effects lead to complex nuclear form 
factors F e m (Q), F S (Q) of Eq. 11091 which makes the relative phase of the 
Primakoff- nuclear interference cross section angle- dependent. Faldt also uti- 
lized a more modern non-uniform nuclear matter distribution. All of these 
theoretical advances were utilized by the Cornell group and their result for 
the 7T° -»■ 77 decay width is T = 7.92 ± 0.42 eV(5.4%)[32]. Taking the 
branching ratio of 1.2% of the Dalitz decay ir° — > e + e~7 into account yields 
r(n°) = 8.21(0.43) • 10~ 17 sec. The error in the Cornell experiment results 
from combining the spread in values obtained using several different kine- 
matic conditions with the systematic uncertainty. Contributions to the latter 
are estimated for the uncertainty in the nuclear shape parameters, the outgo- 
ing pion-nucleon cross sections, accelerator energy, beam luminosity and for 
the maximum opening angle cut. At this level of accuracy, it is worth noting 
that it was assumed that the incoherent cross section is independent of 0(n°). 
This is contrary to the previous assumption, which reduced this contribution 
at small angles due to Pauli blocking. Since its magnitude is determined at 
larger angles where it becomes relatively more important, this might lead to 
a larger contribution under the small angle Primakoff peak and could in turn 
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Figure 7: Left figure: the energy independent Primakoff function f(t) (see 
Eq JllOl and text). Right figure: the center of mass angle in degrees for which 
the Primakoff cross section is a maximum versus photon energy. 

lead to a small reduction of the reported value of r(7r° — > 77). It should 
also be pointed out, however, that the companion 77 measurement [TS] devi- 
ates strongly from the average value of other experiments. Taking these last 
factors into account suggests that the quoted error of 5.4% is possibly under- 
estimated. However, in our opinion, this is the first modern measurement of 
the 7T° lifetime and should be included in a PDG average. 



4.2 Appendix: tt° — > 77 Opening Angle Distributions 

In its rest frame, the 7r° decays into two gamma rays of energy m 7r o/2 with 
an isotropic angular distribution. In the laboratory frame, this distribution 
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Figure 8: Left figure: The Primakoff and nuclear coherent cross sections for 
C for 7r0 production for photon energies of 2 and 5 GeV. Right figure: The 
coherent cross sections (Primakoff, nuclear coherent, and interference) for ir° 
production divided by Z 2 for C and Pb at a photon energy of 5.2 GeV. For 
these figures the parameters for nuclear coherent production are taken from 
the Cornell experiment [32] . 
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Figure 9: Left figure: the laboratory opening angle for the ir° — > 77 . 
Right figure: Recent results for the tt° opening angle [79] from the PrimEx 
experiment [27] . 

is boosted forward towards p^o, with the angular distribution 

dN cos(# 12 /2) 
dfa * sin 2 (fl 12 /2)[sin 2 (# 12 /2) - sin 2 (# 12 _ mm /2) 
sin(6'i 2 _ min /2) = m^/E^o (111) 

where 9 12 is the opening angle between the two gamma rays and O^-min is 
its minimum value. It can be seen that there is a very sharply peaked distri- 
bution in opening angle starting at 6>i 2 _ min . (The singularity at that value is 
removed when the angular resolution of the detector is taken into account.) 
The values of 6 , 12 _ min versus photon energy are shown in Fig. |9l and it can be 
seen that a measurement of the opening angle was able to provide the early 
Primakoff method experimenters with a method to approximately determine 
the photon energy even though they were employing bremsstrahlung beams. 
The opening angle distribution [79] from the recent PrimEx experiment |27] 
is also shown in Fig. |9j This was initiated by tagged photons between 4.9 
and 5.5 GeV, so that Byi-min was between 2.8 and 3.2 degrees. The sharply 
peaked nature of the opening angle distribution is apparent from this figure. 
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4.3 The Direct Measurement 



The most precisely determined lifetime measurement reported in the particle 
data book was performed at CERN in 1985(23]. This was a direct mea- 
surement of the 7T° decay distance at higher energies than the original 1963 
direct measurement discussed in Sec J2[21j . In the 1985 experiment neutral 
pions were produced by 450 GeV/c protons from the CERN SPS incident on 
a tungsten foil[23]. The measurement consisted of two parts: the first was 
a precise determination of the mean decay length; the second part was the 
estimate of the 7r° momentum spectrum. 

To measure the mean decay length, a production target of 70 \i W was 
placed in the 450 GeV/c proton beam. A second foil, placed at a variable 
distance d from the first, was used to convert some of the gamma rays from 
the 7T° decay into electron, positron pairs. Downstream a magnetic spec- 
trometer system detected 150 GeV positrons, and the positron yield Y(d) 
was measured. To illustrate the sensitivity of this experiment to the 7r° life- 
time we calculated Y(d) for mean momentum < p n o >~ 235 GeV[23j. At 
this momentum the relativistic boost factor 7 ~ 1700 and the mean decay 
distance predicted by the axial anomaly is 45/i. The result is shown in FigJTOl 
with the relatively small sensitivity to the value of the tt° lifetime indicated. 
This shows that the experiment had to be performed with excellent accuracy. 
The reported experimental result result was the ratio 

R = [y (250/x) - y(45/i)]/[y(250/i) - y(0)] = 0.3787 ± 0.0078(2.1%). (112) 

The mean decay length A = 46.5 ± 1.0/x(2.1%) was obtained [22] ■ 

To infer the ir° momentum spectrum N(7r°), measurements were made of 
the charged pion momentum distributions 7V(7r =l= ) . In terms of these results 
the 7T° momentum distribution was assumed to be 

N(n°) = kN(tt + ) + (1 -k)JV(tt-). (113) 

Note that only the relative magnitudes of the momentum distributions are 
relevant for the ir° lifetime determination. For the lifetime and systematic 
error determination k was taken to lie in the range 0.50± 0.25 and to be 
momentum- independent. The momentum distributions from that experiment |23] 
are shown in Fig. [TTJ It can be seen that the measurements that were made 
at CERN were all below = 300 GeV/c, the upper limit of the magnetic 
spectrometer |80j. For higher momenta other experiments and estimates were 
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used [231 ED]- The CERN direct experiment utilized forward-produced neutral 
pions (transverse momentum p t ~ 0). Due to a paucity of pion production 
data, the extrapolations utilized experiments performed at larger values of 
the transverse momenta for which the contributions of excited nucleon res- 
onances did not contribute. These resonances contributed to the breaks in 
the estimated value of N(ir + ) shown in Fig{TTJ The extrapolations also uti- 
lized experiments performed on different targets and energies using Feynman 
scaling; they were checked as carefully as possible by comparison with the 
CERN data [80]. 

From the pion momentum distribution N(p n o) the production probability 
of 150 GeV/c positrons N e +(p n o) must be obtained. 

W e +Gv>) = f(Pe,P^)N( P7T o) (114) 

where f(p e ,Pn ) is the probability that the gamma rays from the ir° — > 77 
decay of a pion with momentum p n o produce a positron of momentum p e =150 
GeV/c. This function is evaluated from an integration over intermediate 
photon momenta and is given in Ref. [23] • To a good approximation the 
results for the 7r° lifetime depend only on the average ir° momentum 

j dp T o p n o N e +(p n o) 

< Pn° >= —ri 7 x— (115) 

J dp^o iV e +Qv>) 

With this assumption the mean decay length of the tt° is 

A = < P *° > CT(7r0) (116) 
m w o 

In the evaluation of the experimental results the full momentum distribution 
was used[23]. The use of the average momentum < p n o >~ 235 GeV/c only 
reduced the final lifetime result by 0.8%. 

To obtain a physical understanding of the pion momentum- dependent 
quantities, we made a graph of the individual ingredients in Eqs. 11141 and 
11151 The results are plotted in Fig. [121 which illustrates the sensitivity to 
the high momentum components of the produced pion spectrum iV(7r ). The 
reason is that f{p e ,Pir°) increases with p n o even though N(tt ) is falling. The 
product N e +(p n o) peaks near < p^o >~ 235 GeV/c. 

The quoted experimental result is r(7r°) = (0.892 ± 0.022 ± 0.017) • 10~ 16 
sec, where the first error(2.5%) is statistical and the second(1.9%) is sys- 
tematic [23J. Adding them in quadrature gives a total uncertainty of 3.1%. 
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Knowing the lifetime, the width r(7r°) = Ti/t(h q ) can be obtained, yielding 
r(vr° -4 77) = (7.25 ± 0.22 ± 0.18)eV@. This result is 6.1% (2 a) below the 
prediction of the axial anomaly and 10.3% ( 3a) below the prediction of the 
axial anomaly with chiral corrections. 

In this experiment the momentum distribution of the neutral pions was 
estimated to be the average of the ir + and ti~ momentum distributions as 
indicated by lower energy experiments. As discussed above, the quoted sys- 
tematic error takes the variation in the relative weighting of these two cross 
sections (see Eq. 11131) . and their variation in shape, into account [23] . How- 
ever, one cannot be sure of the accuracy of these assumptions until the 7r° 
momentum spectrum is explicitly measured. 

It is of interest to discuss how seriously to take the discrepancy between 
the direct measurement and the theoreticcal predictions, particularly those 
involving the chiral corrections to the axial anomaly. The measurement 
of the decay distance is very strong, so any problems most likely are in 
the 7T° momentum distribution which was not directly measured. Based 
on Eq. 11161 agreement with the predictions of the axial anomaly (the ax- 
ial anomaly plus chiral corrections) would require an increase of < p n o > 
of 6.1% (10.3%). These seem rather large effects in view of the effort put 
into their determination [23| 180]. Nevertheless, considering how fundamental 
the prediction of r(7r° — > 77) is, a measurement of N(ir°) would be valu- 
able. Fortunately the Compass experiment at CERN is looking into this 
possibility [8Tj. 

4.4 e + e Colliding Beam Measurements 

Another technique uses the e + e~ — > e + e~7*7* — > e + e~P —> e + e~77 reaction 
with colliding beams where P = 7r, 77, 77 , and 7* are almost real photons 
because the final state leptons scatter at small angles and are not detected. 
An experiment was performed at the DORIS II ring at DESY[23] which 
detected the P — > 77 decays in a crystal ball gamma ray detector which 
covered 93% of Air solid angle. Cuts were made on the invariant mass of the 
two gamma rays at the 7r, 77, 77 masses, and the luminosity was measured via 
elastic e + e~ scattering. The efficiency was determined primarily by Monte 
Carlo calculations. The backgrounds due to beam, residual gas interactions 

2 For this article we shall quote the value of r(7r° — s- 77) = BR(ir° —> 77)r(7r°), where 
BR(tt° ->■ 77) = 0.98798(0.032) [TJ. 
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direct : yield versus distance 




Figure 10: Yield versus the separation of the two plates in the direct experi- 
ment at CERN[23]. The three curves are for the values of r(7r° — > 77) which 
correspond to the result of the CERN experiment [23J and the predictions of 
the axial anomaly and chiral calculations as shown in Fig{TJ This figure has 
a supressed zero since ~ 90% of the yield comes from direct production of 
the 150 GeV positrons and is independent of the plate separation d. 
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Figure 11: The momentum distributions of the direct experiment at 
CERN(Fig. 2 of their paper) [23] . These are the spectra produced at 0° by 450 
GeV protons on Ta. The solid curves are the spectra used in the analysis, (a) 
7r + production spectrum: solid circles indicate measured points. (b)7r~ pro- 
duction spectra; solid circles indicate measured points, (c) ir° production spec- 
tra for k = 0.50(see Eq J113p .(d) ir° spectrum which gives 150 GeV positrons. 
The dot-dashed curve and the dash curves represent variations used of the 
7r + spectrum used to estimate the systematic error. Nucleon isobars give 
peaks at the indicated momenta. 
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Figure 12: The distributions of the direct experiment at CERN[23j. The 
(black)dashed curve is the 7r° distribution(normalized to unity at p^o =150 
GeV/c), the (red dashed) curve is f{p e ,V-K°) (which tends to unity as p n o — > 
450 GeV), the (blue) solid curve is p w o N Pe (normalized to unity at it's peak) . 
See text for discussion. 
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and cosmic ray events were eliminated by separate measurements and also 
by stringent cuts on the transverse momenta of the produced meson, S|p t | < 



This method is a generalization of the Primakoff method in which the 
cross section for the 7 + 7* — > P transition is measured. It is unique in that 
it is based on purely electromagnetic physics, whereas the Primakoff effect 
involves a target nucleus. Furthermore it is the only experiment which, when 
carried out at sufficiently high energies, can measure all the widths of the 
pseudoscalar mesons simultaneously. In this case both of the incident gamma 
rays are slightly off shell (i.e., the four momentum transfer q 2 < 0), whereas 
for the Primakoff effect only one photon is slightly off shell. 

For e + e~ collisions the production cross section calculated from QED is 



^p(qlqi) = T(P->77) 16^(0?, + Q2? ~ m 2 P )^F 2 (ql q 2 2 ) (117) 



where mp is the mass of the produced pseudoscalar, q is the three momentum 
transfer of either of the two virtual photons, and F(qf, q%) is the form factor 
for the 7* (of) + 7* (of) — P vertex which is not specified by QED. To make 
estimates of how much F deviates from unity the vector dominance form 

-^(oi > ?1) = (1 ~~ 9i/ m p) -1 )(l ~ 9l/ m p) -1 )) where m p is the mass of the p 
meson, was used. It was estimated by Monte Carlo calculations that the 
stringent cut on H\pt\ restricts < — q 2 >= lOMeV 2 for ir° production (it is 
larger for 77, 77 production) so that the form factor deviation from unity is 
negligible. 

The result for the tt° 77 decay width is T = (7.7 ± 0.5 ± 0.5)eVj2I|- 
This determination does not give the common impression that the error is 
understated. It is important to have a modern version of this purely electro- 
magnetic determination of r(7r° — > 77) with ~ 1% errors. In this connection 
it is also of great interest to perform measurements of the 77 and rj lifetimes 
with comparable accuracy. Fortunately there are groups at Frascati[82j and 
BES[83] looking into this possibility. 




4.5 Appendix: Data Averaging 

The current particle data book average for t(tt ) = (8.4 ± 0.4(5%)) • 10~ 17 
sec. This was obtained from the weighted average of a data set consisting of 
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six experiments which are indicated in FigJTJ The weighted average is 

N N N 

f = ]rr>,/5>,±0>> 4 )- 1/2 

i=l i=l i=l 

Wi = 1/of 

°i = a stat-i + a %s-i ( 118 ) 

It is assumed that there are no correlations between the systematic and 
statistical errors so that they can be added in quadrature. This formula is 
strictly valid for the situation in which there are repeated measurements of 
the same quantity and when the errors are statistical. Consider the case 
where N measurements have been made and all of the total errors = a . 
In this case the resulting error = ao/y/N. As N increases the error can be 
reduced to close to zero! If there are significant systematic errors this is not 
appropriate. 

An example of this is the t] two photon radiative width which is discussed 
in Sec. El In this case there are eight measurements, all using e + e~ collisions 
with total uncertainties which range from 7 to 27%. The results are in 
reasonable agreement with a x 2 /DOF ~ 1. The overall Particle Data Book 
average is T(r] — > 77) = 4.28 ± 0.19(4.4%). This is an example of how 
the calculated error can be arbitrarily reduced due to the fact that there 
are a large number of experiments. The most accurate experiment gives 
r(r/' -> 77) = 4.17±0.10(2.4%)±0.27(6.5%) where the first error is statistical 
and the second is systematic (the dominant error). The weighted average 
is close to this most accurate value, but has a smaller quoted error. It is 
difficult to see how the error can be significantly smaller than that of the 
most accurate experiment. 

In the case of the ir° lifetime there exist significant differences between 
the various experiments. The Particle Data Group handles this situation 
by taking the value S 2 = x 2 /DOF = Eu^r* - T)/(N - 1). If S > 1 
then they multiply the error by S. An example of the problem of using this 
approach consider the history of the results for the 7r° lifetime in the Particle 
Data Book. In the 2008 version [84]. using Eq jl 181 the results for t(tt°) = 
(8.4 ± 0.20) 10 sec. In this case the value of x 2 — 27 with four data points, 
making S 2 = 9.0 which means that there is a very small probability that 
this is correct. In this case the error given by the particle data book is 
multiplied by 3.0 so that the quoted average is t(tt°) = (8.4 ± 0.60) 10 -17 
sec. |84] . In the present on line edition of the particle data book|26j has 
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added two measurements, both of these have larger errors than the 2008 
Particle Data Book average. One is the 1988 e + e~ measurement with a 9 % 
[check] error |24j. The second added measurement is the it — (3 result with a 
13% error [25]. Both of these have a central value which is consistent with 
the previous average so that the x 2 value has not increased. Furthermore 
since two measurements with relatively large errors (small weights) have been 
added the sum. However since two more degrees of freedom have been added 
the sum Sj=i w i is n °t appreciably changed. This means that if S > 1 
the error gets reduced by the square root of N-l, so the new particle data 
book average is now t{ji°) = (8.4 ± 0.4) 10~ 17 sec.[26j. Thus there is a 50 % 
reduction in the error by the addition of two measurements with relatively 
large errors! Clearly this does not seem satisfactory. 

Another point that could be considered is to eliminate any points that are 
far from the average. For t(tt°) this is the 1970 DESY Primakoff point [28] 
which contributes 21.5 out of the 27.0 to \ 2 ■ If this is dropped the average 
value of t(7t°) = (8.7 ± 0.2)10 -17 sec. using the Particle Data Book method. 
Note that this error of 3% is quite small and is comparable to the two best 
measurements [2"3"j |2"T] . Note that these latter two deviate from each other by 
7.5% which is 1.8 times their combined sigmas of 4.2 %. 

The particle data group procedure is a useful way to handle large quan- 
tities of data and get reasonable results. In the case of a quantity for which 
there is a fundamental theoretical interest it is not sufficiently accurate. It is 
clear that some systematic errors have been underquoted, but it is not so clear 
to identify which one. There is no accepted method to combine experiments 
for which the systematic error is dominant or plays a major role. Therefore 
In our view it would be better not to quote an average (particularly for the 
error) under these circumstances. In his case theoretical calculations would 
have to be compared to all of the relevant experiments, not an average value. 
Since one of the main points of this review is that further measurements at 
the ~ 1% level are needed to test this important prediction of QCD, the 
experiments need to be improved. This accuracy cannot be achieved by a 
averaging of experiments with larger errors. 

5 PrimEx Experiment 

After a gap of over two decades there was a new, very accurate measurement 
of the 7T° lifetime, performed using the Primakoff effect, which was recently 
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published [27]. This experiment was performed by the PrimEx collaboration 
in Hall B of Jefferson Lab with 4.9-5.5 GeV photons and was the first Pri- 
makoff effect measurement carried out with tagged photons. The energy was 
determined to ~ 0.1% accuracy for each event. The experiment also ben- 
efited from many technical improvements compared to previous Primakoff 
effect experiments, including a 100% duty cycle (CW) accelerator with a 
smaller beam emittance (a feature that is particularly important for small 
angle experiments), together with greatly improved detectors and electronics. 
The Primakoff effect itself and the previous experiments have been described 
in Subsection 14.11 so that only the specific improvements for the PrimEx 
experiment will be presented below. 

In addition to experimental advances, the improvements over previous 
Primakoff measurements also include advances in the theoretical tools used 
to extract the results. Specifically, the PrimEx experiment utilized recent 
calculations of the nuclear coherent amplitude |76] which represent the strong 
amplitude T s of Eq. 11091 as 

T s = (h- q t )$(0)(F 8 (Q) - wFj(Q)) 
with h = kxe/k, Q 2 = q 2 t + Q 2 min 

q t = p^o sin^o), Q min ~ mlo/k (119) 

— * 

where k is the incident photon momentum, e its polarization vector, and 
$(0) 7^ the forward scattering non spin flip amplitude for photo-pion pro- 
duction on the nucleon. The transverse momentum transfer q t insures that 
the coherent cross section will have the required sin 2 9 dependence. F S (Q) 
is the strong form factor, which takes the final state pion interaction into 
account using the Glauber approach and includes the Faldt correct ion [75], 
which describes the rescattering of the outgoing pions as well as their absorp- 
tion. For light nuclei such as 12 C the next order in the multiple scattering 
series was also taken into account. Such effects were not taken into account 
in the first Primakoff experiments |47[ [28| [33] . The second term in paren- 
theses in Eq J119l accounts for photon shadowing in the initial state. This is 
a two-step process in which the incoming photon produces a vector meson 
(primarily the p ) which in turn produces the emerging 7r°. Fi is the associ- 
ated form factor taking the final state interaction into account The photon 
shadowing parameter w lies between and 1 (none to full shadowing). Fol- 
lowing the empirical evidence the PrimEx analysis assumed w = 0.25 [see 
[75] for references and details]. As part of the systematic error estimate a 
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range of w from to 0.5 was assumed |27j. For the incoherent cross section 
two recent calculations |85[ [HSj were utilized in the PrimEx analysis. One 
used a multiple scattering approach using Glauber theory [85]. The second 
used a cascade model approach [86]. The nuclear incoherent cross section is 
small at small angles (the Primakoff region) and only plays a minor role in 
the extraction of r(7r° — > 77) from the data. 

In a recent theoretical development, photoproduction of 7r°,77,77 mesons 
has been developed in a unified field theoretical basis using photon and Regge 
exchange |8 7]. Full vector dominance has been included. This has been done 
for the Primakoff and nuclear coherent (but not the incoherent) cross sections. 
It is impressive that there is good agreement between the results of this 
calculation, without any empirical parameters, and the PrimEx data and also 
with the fits to the data. In addition, there have been two other calculations 
of the coherent production of the pseudoscalars 7r°, 77, 77 from the proton [88] 
[89] . The second article (89] shows an extensive fit of the Regge parameters to 
the existing data base, which is primarily nuclear production. Both articles 
point out the possibility of performing future Primakoff effect measurements 
with a proton target. This is particularly attractive for the 77 and 77 . 

In the PrimEx experiment incident photons interact with the Coulomb 
field of a nucleus to produce ir° mesons, which quickly decay into two photons 
and are detected in a forward calorimeter as shown schematically in Fig. [T3l 
Tagged photons were used to measure the absolute cross section of small an- 
gle 7T° photoproduction from C and Pb nuclei which make up the target. The 
invariant mass, energy, and angle of the pions were reconstructed by detect- 
ing the decay photons from the ir° — > 77 reaction in a forward calorimeter 
(HYCAL), which was constructed for this experiment. A photograph of this 
detector is shown in FigJTSl The 1152 PbW04 crystals which formed the 
heart of the detector are 2 cm by 2 cm by 20 radiation lengths, which are 
7.5 meters from the target; the published results primarily came from these 
detectors. The outer crystals in HYCAL are Pb-glass. Energy resolution of 
~ 1.5% and angular resolution of ~ 0.02° were achieved. There is an aper- 
ture of 2 by 2 crystals for the highly collimated photon beam. The schematic 
diagram (Fig. [13]) does not show the shielding and the beam dump. These 
and the clean electron beam were sufficient to allow for a low background 
experiment with a sensitive forward calorimeter. The function of the He bag 
was to reduce the number of photons interacting between the target and the 
detector. A plastic scintillator was placed in front of the HYCAL colorime- 
ter to veto charged particles. The magnet, which was placed directly behind 
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the targets, swept the produced electrons away from the detector. Detectors 
were placed behind the magnet to monitor the luminosity using pair pro- 
duction. By turning the magnet to lower fields, measurements of e + e~ pair 
production were made using the plastic scintillator in the coincidence mode. 
By turning the sweeping magnet off, and again using the plastic scintillator 
in coincidence mode, Compton scattering cross sections were also performed. 
To measure the tagging efficiency, HYCAL was moved out of the beam, and 
was replaced by a total absorption counter. This replacement must be per- 
formed at very low currents. The pair production monitors are linear in both 
the higher flux production runs as well as the low flux tagging efficiency runs 
and could be used to interpolate between them. In order to measure the 
7T° production cross sections, precisely measured[90j targets of C and Pb, 
approximately 5% of a radiation length, were used. Further details of the 
experiment can be found in the PrimEx publication [27], on the PrimEx web 
site [35] and in a recent review article by Miskimen|91j. 

To measure the 7 + A — > 7r° + A reaction followed by the ir° — > 7(^1)7(^2) 
decay, where pi,P2 are the four- vectors of the decay photons, two cluster 
events are identified in the HYCAL detector. The center of the hit distri- 
bution is obtained from the distribution of energies in the individual detec- 
tor crystals, and the energy of each photon is determined from the sum 
of these energies. From the inferred energy and angles of the two pho- 
tons the four- vectors Pi,P2 are obtained and the invariant mass distribution 
m 77 = (Pi + V2) 2 = sin 2 (6'i2/2), where ki,k 2 are the energies of two 

gamma rays and 812 is the angle between the two gamma rays, is obtained. 
In addition, one can find the elasticity = (ki + k^)jk^ where k is the inci- 
dent photon energy determined by the tagger to ~ 0.1%. Typical results 
for these quantities in the Primakoff peak region are shown in Fig. HH At 
these small angles the two-photon final states are totally dominated by 7r° 
production as can be seen from the sharp peak for m 77 = m^o ~ 135 MeV. 
The elasticity distribution shows that higher mass states are produced along 
with coherent 7r° production. One important example of such a final state 
is uj meson photoproduction followed by the u — > 7r°7 decay. Another possi- 
ble inelastic mechanism is nuclear excitation or quasi-free meson production. 
The inelastic background is subtracted empirically by fitting the inelastic 
data by empirical polynomial fits. The good angular and energy resolution 
achieved by this apparatus are illustrated in Fig. [14] for the invariant mass 
(~ 1.2%) which allowed identification of the produced ir° mesons with a 
high signal/background ratio. The energy of the emerging pions is also mea- 



59 



sured with good resolution(~ 1.5%) in HYCAL. This resolution (~ 75 MeV) 
does not allow an experimental determination that the pion production is 
coherent, since it allows for the possibility of nuclear excitation. However, 
at small angles this coherent production has been estimated to be small [76]. 
In addition any residual nuclear excitation has been empirically subtracted 
by extrapolating the inelastic background from the measured background at 
higher inelasticities, as shown in FigIT4"l 

The cross sections for forward angle 7r° photoproduction were measured 
on C and Pb targets with 4.9-5.5 GeV tagged photons (having an average 
energy of 5.2 GeV). The resulting experimental cross sections for C and Pb 
are shown in Fig. [TSJ The data were fit by varying the magnitude of each 
of the four contributions of Eq. 11091 — Primakoff, strong, interference, and 
incoherent cross sections. This was done by varying the four parameters 
r 77 . Cg, Cw. & [see [27] for their values]. The resulting cross section fits are 
also shown in Fig. [151 It can be seen that the large forward peak ~ 0.02 deg. 
is dominated by the Primakoff effect, which allows the value of T 77 to be 
accurately extracted. An important test of the dominance of the Primakoff 
mechanism in the small angle region is that the magnitude of this peak scales 
with the nuclear charge as Z 2 . Both the predicted position of the Primakoff 
peak and its separation from the strong 7r° production peak are essential in 
the interpretation of the data. The fact that the theoretical cross sections are 
in such good agreement with the data provides confidence that this separation 
has been done accurately. 

The strong (nuclear coherent) peak would scale in cross section as A 2 (A 
= atomic number) if there were no final state interaction and as A 2 / 3 when 
the mean free path of the outgoing pion is significantly smaller then the 
nuclear radius. In the PrimEx experiment it scales closer to the latter case 
(~ A ' 9 ) and makes the relative magnitudes of of the strong relative to the 
Primakoff peaks smaller in Pb than in C. The angle for which the the strong 
cross section peaks is smaller in Pb than in C, due to its larger radius (which 
increases ~ A 1 / 3 ). The strong-Primakoff interference cross section is the only 
nuclear contribution near the Primakoff peak (~ 0.02°) and enters at the few 
% level. Its strength reflects the positions and magnitudes of the strong and 
Primakoff cross sections. Thus the values of the radiative width T 77 obtained 
from the C and Pb data pose a stringent test on the model-dependence of 
the result. Consistent results for T 77 were obtained from the data for these 
two targets, supporting the idea that the small (few %) nuclear effects being 
subtracted under the Primakoff peak are well described by the theoretical 
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Figure 13: Schematic layout of the PrimEx experimental setup showing 
the incident electron beam, photon energy tagging system, target, sweeping 
magnet, He bag, electron pair spectrometer, veto counter, and HYCAL de- 
tector; shown in more detail in the insert. It consists of an inner section of 
PbW04 crystals and an outer section of Pb-glass detectors. 
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Figure 14: Observed invariant mass(left) and inelasticity distribution data. 
See text for discussion. 

calculations whose magnitudes are fit to the larger angle data. 

The largest two systematic errors of the experiment are 1.0% for the 
photon flux and 1.6% for the yield extraction. The uncertainty in the flux 
is caused by instabilities in the photon beam and detection. The error in 
the yield extraction is primarily due to uncertainties in the background sub- 
traction, as illustrated in Fig. [TH The total error of 2.1% was obtained 
by summing all of the errors in quadrature, since there are no known cor- 
relations between them. As stated above, this is the smallest systematic 
error for a photon-induced cross section that we are aware of. The ability 
to accurately measure cross sections with this apparatus was tested by mea- 
surements of the Compton effect and pair production, which are accurately 
predicted. The measurements agree with theory to < 2%, which is better 
than the systematic error for the Primakoff effect. 

The PrimEx data were independently analyzed by several groups in the 
collaboration. This included the event selection and independent develop- 
ment of software, apart from sharing the same photon flux routines which 
were independently checked by measurements of the pair production and 
Compton cross sections. The results for the 12 C and 208 Pb targets are shown 
in Fig. [T6j It is gratifying that these two independent analyses agree with 



62 



9(n), (degree) 



9(h), (degree) 



Figure 15: Cross sections da/d8 in /xbarn/rad versus the lab pion angle for C 
(left)and Pb(right) at an average photon energy of 5.2 GeV. The individual 
contributions were obtained by a fit to the data (see text for discussion). 
The Primakoff contribution peaks ~ 0.02° , the strong contribution (red dot- 
ted curve) peaks ~ 1.6° in C and ~ 0.8° in Pb with a smaller secondary 
maximum~ 1.8°. The interference contribution (green dotted curve) peaks 
~ 0.9°, 0.3° in C and Pb. The incoherent contribution(solid blue curve) rises 
gently with 9 n o. 
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each other within the errors. Perhaps even more important is the fact that 
the width extracted from the C and Pb targets agree within errors! This 
verifies the Z 2 dependence of the Primakoff cross section at the few % level 
and is a strong indication that all of the nuclear effects have been properly 
taken into account. The PrimEx result is r(7r° -»■ 77) = 7.82 ± 1.8% ± 2.2% 
where the first error is statistical and the second is systematic. Combining 
them in quadrature gives a total error of 2.8%. This result for T 77 is within 
one standard deviation of the theoretical prediction [29 [ |30| |3T] and most of 
the results of previous measurements [14"] as shown in Fig. [TJ 

The Primakoff experiment represents an effort to reduce the experimental 
error in r(7r° — > 77) to the 1% level which is the present theoretical accuracy. 
Because of the this requirement the PrimEx collaboration proposed and ex- 
ecuted a follow up experiment in the Fall of 2010 at JLab. The goal is to 
reduce the total error to less than 2%. Targets of 28 Si and 12 C were employed 
and the statistics were significantly improved relative to the original PrimEx 
run. The data analysis is presently in progress. 

We conclude this section with a more up to date summary of the present 
status of the r(7r° — >• 77) decay width is presented in Fig. [T7] and Table 
|2j Following the recommendation that was made in Sec. H] several of the 
previous Primakoff measurements which were performed at lower energies, 
and analyzed with imprecise theoretical models are not included. Following 
the suggestion made in Sec. !4.5l no average value is given, just the individual 
data points. As we reach a new level of precision for the 7r° lifetime measure- 
ment the situation is not completely satisfactory. The two most accurate 
experiments are not in agreement. The direct experiment |23j differs from 
the PrimEx result [27] by 7.5% which is 1.8 times their combined sigmas (in 
quadrature) of 4.2 %. This clearly needs further investigation. 

In summary three of the four experiments are in agreement with both the 
axial anomaly and the chiral predictions. They are clearly not of sufficient 
accuracy to show the 4.5% increase predicted in the width by the chiral 
corrections and to fully test them. 

6 Related Measurements 

There are two related extensions of the physics of the ir° — > 77 rate. One is 
to consider the decay rates of the rj and 77 pseudoscalar mesons. The other is 
to consider when one or two of the photons are off shell, ie at a finite value of 
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Figure 16: The r(7r° — > 77) widths measured by the PrimEx experiment [27] 
for the C and Pb targets, along with the chiral predictions |29l 1301 I3"T] . 
The errors are statistical and systemic added in quadrature. The lower hor- 
izontal line is the predictions of the chiral anomaly (31 H] (r(vr° — > 77) = 
7.760eV, r(7r°) = 0.838 • 10~ 16 sec). The upper horizontal lines are the chiral 
predictions^ EES EI](r(7r 77) = 8.1eV,r(7r°) = 0.80 ■ 10" 16 sec.) with 
it's estimated 1% error. 
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Figure 17: tt° — > 77 decay width in eV (left scale) and t(7t°), the mean 7r° 
lifetime in units of 10~ 16 sec. (right scale) measured by the most accurate 
experiments. These include the 1974 Cornell Primakoff measurement [32] . 
the 1985 direct measurement at CERN[23j. the 1988 e + e~ experiment at 
DESYpM] . the 2011 PrimEx experiment at JLab[27]. The lower dashed line 
is the predictions of the chiral anomaly[HJ|I] (r(vr° — > 77) = 7.760eV, t(ti ) = 
0.838- 10~ 16 sec). The upper solid line is the chiral prediction and the dotted 
lines show the estimated 1% error [291 E3 EI] (r(7r° 77) = 8.10eV, t(tt°) = 
0.80 ■ 10 -16 sec). For the relationship between r(7r° — > 77) and r(7r°) see 
Eqffl 
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Table 2: Most Accurate Measurements of r(7r° — » 77) and t(tt°). For their 
relationship see Eq{TJ 



Reaction 


Reference 


r(vr° -> 77) eV 


r(^°)/10- 16 sec 


Primakoff:1974 


m 


7.93(0.39) 


0.820(0.040) 


direct: 1985 




7.25(0.22) 


0.897(0.028) 


e+e-:1988 


m 


7.70(0.71) 


0.845(0.078) 


Primakoff:2011 




7.82(0.22) 


0.832(0.023) 



the four momentum transfer carried by one of the two photons. This could 
be either space like (q 2 < 0) or time like(g 2 > 0). The space like transitions 
are accessed by the e + e~ —> e + e"7*7* — > e + e~P reactions (P = 71,77,77) 
when one of the leptons in the final state are detected at non-zero angles. 
This experiment has been carried out at DESY for ir°, rj, rj production by the 
CELLO collaboration for q 2 values from approximately -0.3 to -3 GeV 2 [92\. 
The results were compared to the dipole form factor F(q 2 ) = 1/(1 — q 2 /A 2 ) 
where A is fit to the data. To first approximation the three form factors 
could be fit with A ~ m p (the p meson mass), although for the best fits 
A did differ for each meson. These differences can be explained by ChPT 
calculations [93]. The results for the transition radius of the 7r° is 0.65± 0.03 
fm [92], close to the RMS radius of the of the charged pion. These results 
for the transition radius are model dependent since they are performed at 
relatively high values of q 2 . A new measurement for the 7r° transition form 
factor at low q 2 values has been proposed at Frascati[82] with the KLOE-2 
e + e~ colliding beam detector. It is also possible to probe the low q 2 region 
by studying virtual Primakoff meson production eA — > e tt A|94]. 

For the time like region low q 2 measurements have been performed by 
observing the tt° — > 7e + e~ reaction using neutral pions produced in the 
7T~p — )■ 7r°p reaction with stopping pions [95j [96]. By expanding the transition 
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form factor at low q 2 these results are consistent with those measured in 
the space like region, but with larger errors. The radiative corrections to 
this process have been worked out in detail in anticipation of more accurate 
measurements [97] . 

At the present time the measurements of the T(r],r] — > 77) rates have 
significantly larger errors than for the 7r°[Il]. These have come from e + e~ 
experiments. For the rj the Particle Data Book uses four measurements 
carried out between 1985 and 1990. They are all in good agreement and the 
resulting average is r(?7 -» 77) = 0.510 ±0.026(5.1%)keV[H]. For the rj the 
Particle Data Book uses eight measurements carried out between 1988 and 
1998. They are also in good agreement and the resulting fit is T(r] — > 77) = 
4.30 ±0.15(3.5%) keV[H]. 

There has only been one Primakoff measurement for the rj. This has re- 
sulted in a significantly smaller value of r(rj -> 77) = 0.324±0.046(14.2%)keV[78] 
and it was not used in the Particle Data Book fit[T3]. This was performed 
at Cornell with bremsstrahlung beams with end points between 5.8 and 11.5 
GeV with five targets between Be and U. It is of interest that this was 
the same group that performed such an accurate measurement for the ir° 
lifetime [32]. The reason that Primakoff measurements of heavier mesons 
is much more difficult than for the 7r° can be understood on the basis of 
the Primakoff effect discussed in Subsectior j4.ll From Eq. 11091 it can be 
seen that the differential cross section for the Primakoff effect for the eta 
meson is proportional to k 4 Z 2 T(r] — > 77) /mL If we note that the width 
T(r) — > 77) oc w? n then the Primakoff cross section decreases oc 1/m 4 , so 
that it is considerably smaller than for the ir°. In addition the Primakoff 
peak position 9 max = m 2 /(2k 2 ), which is considerably larger than for the 
7T°. The combination of these two effects means that the nuclear interfer- 
ence is considerably more of a problem for the rj than for the 7r°. This can 
be seen by looking at the figures for the yields for 77 production [78]. In 
a recent publication the nuclear incoherent effects have been re-evaluated 
and the conclusion is that the width measured by the Cornell group was 
changed to T(r] — > 77) = 0.476 ± 0.062(13.0%)keV[9"8"]. in reasonable agree- 
ment with the PDB fit[T4] based on e + e~ two photon experiments. Although 
this reanalysis is a significant improvement of the treatment of the incoherent 
nuclear rj production, it was not meant to replace a new experiment. The 
PrimEx collaboration has an approved experiment to measure the Primakoff 
r\ production cross section from the proton with the upgraded 12 GeV JLab 
facility [92], with the goal of reaching 2.2% accuracy. 
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Improving the precision of the 77 and 77' two photon decay rates is im- 
portant for several reasons. It is important in determining the 77,77 mixing 
which enters their lifetime calculations as well as that of the 7r°. It is also 
important since all of the decay rates are based on the well known branch- 
ing ratios and the two gamma decay rates[14J. Determination of the isospin 
breaking 77 — > 3n decay rate provides an independent determination of the 
mass difference of the up and down quarks m^ — m u . Finally there is the more 
speculative issue of the nature of the 77 meson, which has too large a mass to 
be a Nambu-Goldstone Boson, but is so in the large N c limit. The question 
of its glueonic content has been a long standing issue. Finally we note that 
the masses and mixing of the 77, 77 mesons have been recently calculated on 
the lattice [TOO]. 

7 Future Outlook 

The 7T° — > 77 decay is perhaps the best example of a process that proceeds 
primarily via the chiral anomaly [HJ H], which is an essential component of 
QCD. 

This possibility exists due to spontaneous chiral symmetry breaking, 
which makes the ir° the lightest hadron, and consequently its primary decay 
mode is into two gamma rays. The chiral anomaly represents breaking by 
the electromagnetic field of the symmetry associated with the third isospin 
component of the axial current [3j H] . This ir° decay provides the most sen- 
sitive test of this phenomenon of symmetry breaking due to the quantum 
fluctuations of the quark fields in the presence of a gauge field. Consider- 
ing the fundamental nature of the subject, and the 1% accuracy which has 
been reached in the theoretical lifetime prediction, it is important for future 
experiments to aim for a comparable level of precision. 

The prediction of the chiral anomaly for the ir° lifetime has no adjustable 
parameters and is exact in the chiral limit — (m^m^m,,-) — > 0. The chiral 
corrections involve isospin breaking and mix the 77,77 mesons into the 7r°; 
consequently, they are proportional to m u — m^. There exists a (4.5 ± 1.0)% 
increase to the n° decay rate predicted by the chiral calculations [221 EDI EI] 
relative to the anomaly. Although the theory is consistent with experiment 
at the present level of accuracy, there is a significant need for future ex- 
perimental improvements to fully test the chiral prediction. At the present 
time theory is ahead of experiment in that the estimated theoretical error 
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of 1% in the chiral calculations [221 EDI EI] is significantly smaller then the 
experimental uncertainty. The theoretical error arises from uncertainties in 
the low energy constants and is not due to higher orders in the chiral ex- 
pansion. Therefore it is not likely that the present theoretical error can be 
significantly reduced using chiral perturbation theory. On the theoretical 
side further progress probably requires lattice calculations. A first lattice 
calculation of rj, r\ mixing has recently been reported |lUUj . 

It is very important to have modern experiments performed with all three 
techniques at the 1% level. The PrimEx experiment has a quoted accuracy of 
2.8% [27] and the direct experiment at CERN 3.1% . However, the difference 
between the central values of these results is 7.3% . This discrepancy clearly 
needs to be resolved if further progress is to be made. The PrimEx group 
has had a second run at JLab using 12 C and Si targets with improved 
systematics. Using the analysis techniques developed for the experiment 
the goal is to reduce the error to < 2%. We know of plans to remeasure the 
direct experiment with the Compass experiment at CERN[81J. There are also 
plans to remeasure the ir° lifetime at Frascati[S2] and at Belle [83]. This latter 
effort can also measure the rates for the t] and rj as well as the q 2 dependence 
associated with Dalitz decay. We look forward to future developments. 
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